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ABSTRACT 

We present mid-infrared spectra of T Tauri stars in the Taurus star-forming region obtained with 
the Spitzer Infrared Spectrograph (IRS). For the first time, the 5-36 //m spectra of a large sample 
of T Tauri stars belonging to the same star-forming region is studied, revealing details of the mid- 
infrared excess due to dust in circumstellar disks. We analyze common features and differences in the 
mid-IR spectra based on disk structure, dust grain properties, and the presence of companions. Our 
analysis encompasses spectral energy distributions from the optical to the far-infrared, a morphological 
sequence based on the IRS spectra, and spectral indices in IRS wave bands representative of continuum 
emission. By comparing the observed spectra to a grid of accretion disk models, we infer some basic 
disk properties for our sample of T Tauri stars, and find additional evidence for dust settling. 
Subject headings: circumstellar matter — planetary systems: protoplanetary disks — stars: pre-main 
sequence — infrared: stars 
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1. INTRODUCTION 

The Taurus- Auriga star-forming region is a relatively 
nearby (140 pc: lKenvon et alll994bUBertout et al.ll999J) 
molecular cloud complex, which is characterized by lower 
extinction and more isolated, low-mass star formation 
than other star-f orming regions at comparable distances, 
like Ophiuchus ijKenvon fc Hartmannl 119951 hereafter 
KH95). It is therefore fairly well-studied from wave- 
lengths ranging from the X-rays to the radio and used to 
test theoretical predictions of the star formation process. 
Given the young age of the Taurus complex (1-2 Myr), 
many of its pre-main-sequence stars are still surrounded 
by material left over from their formation process and 
accrete gas and dust from circumstellar disks. 

Dust grains in these disks are heated by irradiation 
from the star and from the accretion shocks at the 
stellar surface, as we ll as by local viscous dissipation 
ijD'Alessio et al.lll998j) . and emit excess emission, i.e. at 
flux levels above those expected from a stellar photo- 
sphere, at infrared wavelengths. Most disks have flared 
disk surfaces, i.e. the disk thickness increases more 
steeply than just linearly with radius; this flaring allows 
disks to absorb more radiation from the central source 
tha n a flat, thin disk, and thus to become hotter (see, 
e.g.. lKenvon fc Hartmanr]ll987j) . 

Over time, accretion rates decrease, and, while the disk 
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material is depleted, some of the dust and gas could be 
built into planets. Young, low-mass (< 2M Q ) stars sur- 
rounded by these accretion disks are referred to as classi- 
cal T Tauri stars (CTTSs). Once accretion has virtually 
stopped and disks are either passive, i.e. only repro- 
cessing stellar radiation, or already dissipated, a young 
star is called a weak-lined T Tauri star (WTTS). The 
distinction between CTTSs and WTTSs is based on the 
equivalent width of their Ha emission line (EW[H«]): 
CTTSs have broad and asymmetric Ha line profiles, 
which are generated in the magnetospheric accretion 
flows ijMuzerolle et al J 12003^1 . resulting in an EW[Ha] 
larger than 10 A, while WTTSs have an EW[Ha] of less 
than 10 A; in WTTSs, strong surface magnetic fields pro- 
vide the energy for strong X-ray, UV, and Ha emission. 

Based on an evolutionary sequence derived from the 
shape of the spectral energy distribution (SEP) in the 
infrared lILada fc Wi 1 kin dll9^ 171^119^ . young, low- 
mass stars are divided into three classes. Class I objects 
have rising SEDs over the infrared spectral range, while 
Class II objects have either flat or somewhat decreasing 
SEDs; Class III objects have little infrared excess, and 
their SEDs can be understood as reddened stellar photo- 
spheres for wavelengths < 10 /im. Thus, Class II objects 
are characterized by the presence of mid-infrared excess 
emission due to dust in circumstellar disks, while Class 
III objects have, for the most part, already dissipated 
their disks, and therefore their infrared flux is close to 
photospheric levels. Classical T Tauri stars fall into the 
Class II category; weak-lined T Tauri stars can be con- 
sidered Class II or Class III objects, depending on the 
presence of disks around them. When a young star has 
almost completely dissipated its circumstellar material, 
it is usually classified as a Class III object. Therefore, 
all Class III objects are WTTSs, but not all WTTSs are 
Class III objects, as defined by their infrared SEDs. 

Mid-infrared spectra of T Tauri stars (mostly CTTSs) 
reveal the details of the thermal emission by dust grains 
in the inner regions of circumstellar disks, from a few 
tenths to several AU from the central object. By de- 
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riving the structure and composition of the dust in pro- 
toplanetary disks, we can study the conditions and pro- 
cesses leading to planetary formation. In particular, dust 
growth and settling are seen as key steps in the formation 
of larger bodies. 

The Spitzer Space Telescope l)Werner et al]l2004|) en- 
ables us, for the first time, to obtain spectra from 5 to 
36 /im of T Tauri stars covering a variety of flux den- 
sities, from tens of mJy to tens of Jy. As part of a 
larger Infrared Spectrograph 6 (IRS; iHouck et alJ 12004^ 
guaranteed-time observing program, we observed 150 
young stellar objects in the Taurus- Auriga star-forming 
region, ranging in evolutionary state from the protostel- 
lar, or Class I, stage, when a forming star is still en- 
shrouded by an envelope, to the Class II and III stages, 
when disk material is accreted and dissipated. 

Our Taurus sample is largely based on the sample stud- 
ied by KH95, who compiled and analyzed photometric 
data of known young stellar objects in the Taurus- Auriga 
star-forming region. Here we present 85 IRS spectra of 
Class II objects and 26 spectra of Class III objects; the 
remaining spectra of the Class I objects in our sample 
will be presented in a subsequent paper. We analyze the 
spectral energy distributions and IRS spectra of Class II 
and III objects in terms of common features and pecu- 
liarities arising from their disk structure, the properties 
of dust grains embedded within, and their surrounding 
environment. The large-scale environment has little im- 
pact on objects in Taurus, unlike young stars in OB as- 
sociations or young clusters, where the cluster members 
are subject to the influence of high-energy radiation and 
powerful winds from the OB stars. 

In § we introduce our Spitzer IRS observations and 
data reduction; in §|3we present spectral energy distri- 
butions, from the optical to the far-infrared, of the Class 
II and III objects in our sample, as well as a morphologi- 
cal sequence of the IRS spectra of our 85 Class II objects; 
in § 0] we compare our data with a set of accretion disk 
models to infer some disk properties, present a new me- 
dian of mid-infrared fluxes of Class II objects in Taurus, 
and estimate upper limits for dust masses around our 
Class III objects. We discuss some of our results in § 
and give our conclusions in § 

2. OBSERVATIONS AND DATA REDUCTION 

All but three of our Taurus targets were observed with 
the IRS on Spitzer during IRS observing campaigns 3 
and 4, from 2004 February 6 to 8, and from 2004 Febru- 
ary 27 to March 5. The remaining three targets (FQ 
Tau, GK Tau, and V830 Tau) were observed during IRS 
campaign 12, from 2004 August 30 to 31. For all tar- 
gets the full mid-IR spectrum from 5 to 40 /im was ob- 
tained by either using the two low-resolution IRS mod- 
ules (Short-Low [SL] and Long-Low [LL], 5.2-14 /im and 
14-38 jum, respectively, A/ A A ~ 90) or the SL module 
and the two high-resolution modules (Short-High [SH] 
and Long-High [LH], 10-19 /im and 19-37 /im, respec- 
tively, A/ A A ~ 600). The high-resolution modules were 
usually only used for observations of bright (> a few 
Jy) targets. While some targets were observed in IRS 

6 The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation funded by NASA through the 
Jet Propulsion Laboratory and the Ames Research Center. 



staring mode, most of them were observed in mapping 
mode; with the latter, a 2x3 step map on the target was 
carried out, with three steps separated by three-quarters 
(for SL) or half (for the other modules) of the slit width 
in the dispersion direction and two steps of one third of 
a slit length in the spatial direction. For staring mode 
observations, the target was placed subsequently in two 
nod positions along the spatial direction of the slit, at 
1/3 and 2/3 of the slit length, respectively. 

As a first step in the data reduction process, we fixed 
all bad pixels in the arrays (identified by the pixel masks, 
and including the so-called "rogue pixels" in SH and LH) 
using a simple interpolation over neighboring, good pix- 
els. We then extracted and calibrated our spec tra using 
the SMART software tool <|Higdon et all 12004(1 . When 
taking low-resolution spectra, the target is first placed 
on the slit yielding the 2nd order spectrum and then on 
the slit yielding the 1st order spectrum; each time both 
orders are observed. We subtracted the sky background 
by using observations taken at the same nod position, but 
in different orders. The spectra were extracted using a 
variable-width column extraction, which varies with the 
width of the point-spread function, and ca librated with 
a La c (Al V) and its template spectrum l|Cohen et alJ 
|2003|) . For the high-resolution spectra, the spectra of tar- 
get and calibrator, £ Dra (K2 III), were extracted using 
a full slit extraction. 

Even though we were not able to subtract any sky 
background from our high-resolution spectra, the contri- 
bution from any background emission was only notice- 
able in our LH spectra. To correct for this contribution, 
we applied a scalar correction to each LH spectrum to 
achieve a smooth transition at 19 /im, where SH and LH 
meet. Our SH spectra usually matched the flux level of 
our SL spectra in the overlap region between SL and SH 
(10-14 /im); given the smaller aperture and shorter wave- 
length coverage of SH, any background contribution is 
expected to be negligible. All of our high-resolution data 
presented in this paper were rebinned to a resolution of 
A/AA = 200, and the SH data were truncated below 13 
/im. The rebinning of our SH and LH spectra resulted 
in a clearer representation of the shape of the continuum 
and in a more homogeneous set of spectra, independent 
of the modules used for A > 14 /im. 

For mapping mode observations, we extracted and cal- 
ibrated all map positions, but the final spectrum was 
usually obtained by averaging the two observations in 
the map's center. For staring mode observations, the 
final spectrum is an average of the two nod positions. 

We estimate our absolute spectrophotometric accuracy 
to be of the order of 10%, in some cases as good as 5%. 
The scatter of neighboring flux values determines our 
relative accuracy and thus our ability to reliably iden- 
tify spectral features. This relative accuracy depends 
not only on the brightness of the object, but also on the 
modules used: while in SL, and usually also in LL2 (14- 
21 /im) and in SH, spectral features above the noise level 
can be considered real, the LL1 (20-36 /im) and espe- 
cially the LH spectra are often dominated by artifacts 
due to unresolved calibration issues, which makes the 
identification of real features difficult. Therefore, spec- 
tral features beyond 20 /im need careful checking before 
they are identified. In addition, the wavelength region 
between 13.5 and 14.5 /im, where SL and LL meet, is 
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affected by larger uncertainty in the calibration at the 
order edges and by possible flux mismatches between SL 
and LL. In some cases a spectral feature seems to appear 
in this wavelength region; at this point, we are not sure 
whether this is a real feature or just an artifact intro- 
duced by the order edges. 

It should be noted that many of our targets are mem- 
bers of close multiple systems, ranging in separation from 
tenths of an arcsecond to several arcseconds (see Tables 
Q]and|3I). Given that the narrowest IRS slit width is 3'.'6 
(SL), several targets could not be resolved by our obser- 
vations. For these targets, we show the combined spec- 
trum of the unresolved components; however, in many 
cases one component clearly dominates the mid-IR spec- 
trum. For a small group of targets, a companion entered 
only partially in the wider slits (4'.'7 for SH, 10'.'6 for LL, 
11" 1 for LH) and more fully at the longer wavelength 
ends of these modules. For these objects we made sure 
that the contribution of the companion was small before 
interpreting the object's spectrum, or we excluded the 
"contaminated" part of the spectrum from the analysis. 

A few objects presented problems in their reduction; 
they will be listed below. 

IT Tau. This classical T Tauri star is a 2'.'4 binary. Both 
components entered the SL and LL slits, even though 
some flux of the B component was lost at longer SL wave- 
lengths due to the larger size of the point-spread function. 
The A component should dominate the short-wavelength 
part of the IRS spectrum; however, it is likely that in LL 
(> 14 /urn) A a nd B contribute about eq ually, since IRAC 
measurements IfHartmann et aLll2005a|) indicate that the 
contribution of the B component to the total flux of the 
system increases from 15% at 3.6 /xm to ~ 30% at 8 /xm. 
This is also suggested by the IRS spectrum, where we 
multiplied SL by 1.25 to match it with LL, thus creating 
a combined spectrum for A and B. Due to this scaling 
the SL spectrum of IT Tau carries a large uncertainty. 
CW Tau, V773 Tau. These two objects are part of a 
mispointed mapping-mode observation. We created con- 
tinuous IRS spectra by using the brightest map positions 
for each star; however, due to likely slit losses the abso- 
lute flux levels and, to a lesser degree, the shape of the 
spectra could be off. 

CX Tau, FN Tau, FO Tau. The spectra of these objects 
were also affected by mispointed observations, and there- 
fore their IRS spectra are less accurate. 
MHO 3. This young star is so bright that a small part of 
its LL spectrum (from about 22 to 26 /im) was saturated. 
AB Aur, RY Tau. These objects are also very bright and 
partly saturated the central map position in SL. We re- 
covered an unsaturated spectrum by using an off-center 
map position for SL and scaled it to match the SH and 
LH observations. 

V892 Tau. Due to the brightness of this star, SL and LH 
were partly saturated. Similarly to RY Tau, we used off- 
center map positions for SL and LH and matched them 
to the central map observation of SH. 

3. THE SAMPLE 

3.1. Characterization and Classification 

The targets of our Taurus sample were selected from 
the objects analyzed by KH95, and from objects in the 
IRAS Faint Source Catalogue, selected on the basis of 
their infrared colors. For the latter, we used a filter on 



the brightness and on the 12/25, 25/60, and 60/100 /im 
flux ratios (the first one being the most important one): 
12/25 < 1.0, 25/60 < 1 . 2, an d 0.1 < 60/100 < 1.5, 
similar to Beic hman et al.l (^992) . Sources from the Faint 
Source Catalogue that do not have common names are 
labeled with their IRAS identifier, e.g. F04101+3103. 

Spectral indices have commonly been used to clas- 
sify the SEDs of protostars and pre-main-sequence stars 
ljLadall1987t I Adams et a,l]IT987^ . The slope of the SED 
from the near- to the mid-infrared (2.2-25 /im) in a log- 
log plot is used to define Class I, II, and III objects: 
n = d\og(\F\)/d\og(\), where n > identifies a Class I 
object, — 2 < n < a Class II object, and n < — 2 a Class 
III object. A value of n of —3 is equivalent to a stellar 
photosphere. However, this is just a simple classification 
scheme, since it does not take into account the detailed 
shape of the SED in this wavelength range; it describes 
the shape of the SED from 2 to 25 /im as a simple power 
law XF\ oc A". For example, some objects identified as 
Class I objects might actually be Class II objects with 
cleared-out inner disk regions, which are characterized by 
little excess emission at near-IR wavelengths and a large 
mid-IR excess, resulting in a steep slope of the near- to 
mid-infrared SED. Only mid-IR spectra reveal spectral 
features that can aid in the detailed evolutionary classi- 
fication of a young stellar object, like the presence of the 
10-/im silicate emission feature, which is indicative of a 
flared circumstellar disk. 

We separated our sample into Class II and III sources 
based on the existence of an infrared excess in the IRS 
spectral range. This subdivision resulted in a reclassi- 
fication of some objects previously identified based on 
their SED slope between 2 and 25 ptm. Since the 25- 
/im data points were usually derived from IRAS mea- 
surements with limited spatial resolution and lower sen- 
sitivity, many previous classifications were more uncer- 
tain. Reclassification was done for FP Tau, V836 Tau, 
VY Tau, and ZZ Tau, which were originally classified as 
Class III objects (e.g., KH95), and for DI Tau and HBC 
423, which were considered Class II objects. Also, we 
reclassified a few Class II objects as Class I objects due 
to the presence of ice and silicate absorption features in 
their IRS spectra, as well as their SED shape; among 
these objects are CoKu Tau/1, HL Tau, and LkHa 358, 
which will be included in a subsequent paper. 

Accretion signatures, like broad Ha emission and an 
excess in the UV or blue spectral range, are an indica- 
tion for the presence of inner disk material. Combined 
with magnetospheric accretion models, these signatures 
are u s ed to derive mass accre ti on rates l| G i ullbrin^^lJ 
119981 lHartmann et all Il998t iMuzerolle et all I2003a1. 
Classical T Tauri stars in Taurus have average mass 
accretion rates of ~ 10~ 8 M Q yr -1 ; accretion rates 
seem to increase w ith the mass of the star and also 
decre ase over time llCalvet et alJ 120041: IMuzerolle et alJ 
120051 lHartmann et al.lll998|) . 

However, Ha emission is not necessarily correlated 
with the presence or absence of a disk; several stars 
identified as WTTSs do have substantial disk masses, 
as derived from the ir millimeter emission (see, e.g., 
iBeckwith et al.lll990|) . or have considerable mid-IR ex- 
cesses, as our sample shows. This observational fact 
could be explained with the presence of inner disk holes: 
a star with a missing inner disk might not be accreting, 
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and thus be classified as a WTTS, bu t still possess a 
prom inent outer disk (e.g., CoKu Tau/4: lD'Alessio et alJ 
2005). On the other hand, all stars with inner disks, 
as inferred from their infrared excesses, are accret- 
ing and t hus CTTSs. A few CTTSs (e.g., TW Hya 
GM Aur: I.Tavawardhana et aUllQQflt lUchida et al.ll2004t 
ICalvet et al.l2005|) have low near-IR excesses, but promi- 
nent mid-IR excesses, which could be an indication of 
grain gro wth in the inner dis k. An extreme example is 
DM Tau ijCalvet et all 12005(1 . which is an accreting T 
Tauri star but has no near-IR excess. 

Many, if not most, T Tauri stars in Taurus have com- 
panions. A list of our Class II objects belonging to 
multiple systems can be found in Table ^ many com- 
ponents are separated by fractions of an arcsecond, and 
usually they form a binary or a hierarchical quadruple 
system. The multiplicity of young stars adds complexity 
when their stellar and circumstellar properties are stud- 
ied. Some syst ems are surround ed by circumbinary disks 
or rings (e.g.. IKrist et al.|l2002|) . often in addition to a 
circumprimary and a circu msecondary disk. These disk s 
can be tidally truncated l|Artvmowicz fc Lubowl [19941; 
in general, disks are found well within or outside binary 
orbits. Therefore, the components of a binary system can 
be either both CTTSs or both WTTSs, especially in close 
(< a few arcseconds) binary systems, or form a mixed 
pair llPrato fc SimoullT997t IDuchene et a,l.llT99l. A cir- 
cumbinary reservoir can replenish both the circumpri- 
mary and circumsecondary disk; for wider companions, 
the circumprimary disk is preferentially replenished, re- 
sulting i n a stronger infrared e xcess for the primary com- 
ponent i White fc GhedfeOOlj) . Depending on the pres- 
ence of material in its inner disk and thus its accre- 
tion signatures, the primary may be classified as CTTS 
or WTTS; eventually, both components will evolve into 
Class III objects (and thus WTTSs), once all the disk 
material has dispersed. 

3.2. Spectral Energy Distributions 

Our full sample of Class II and III objects is listed in 
Tables and 01 respectively, together with the adopted 
spectral type, visual extinction Ay, and T Tauri star 
type ( CTTS, WTTS) taken from the literature (mainly 
KH95: iWhite fc Ghedl200lD . Due to the variability of 
these sources and their multiplicity (where applicable) 
some measurements are quite uncertain. Several spec- 
tral types are not known or are only determined with an 
uncertainty of a few subclasses, and ages and mass ac- 
cretion rates depend on the models used to derive them, 
which can cause uncertainties of an order of magnitude. 
When derived using the same assumptions, these quan- 
tities usually agree within a factor of two. When sepa- 
rating the star s of our sample into C TTSs and WTTSs, 
we adopted the lWhite fc Basril pOO^ criteria for CTTSs, 
where the traditional EW[Ha] boundary of 10 A between 
WTTSs and CTTSs is adjusted for spectral type, with 
a delimiting width smaller than 10 A for spectral types 
earlier than K7 and a larger width for spectral types later 
than M2.5. 

For each object with known spectral type, we checked 
whether the Ay value from the literature or the Ay value 
we derived from observed V-I colors yielded dereddened 
optical photometry consistent with the intrinsic colors of 
a main-sequence star with the same spectral type as the 



object under consideration. If the dereddened flux values 
were off, we adjusted the Ay values, usually by less than 
a factor of two, until a value was found that brought the 
observed and expected colors into agreement. Several 
Ay values remain uncertain, but since for most Taurus 
targets the extinction is low, changes in Ay by less than 
a factor of two will only minimally affect the shape of 
the IRS spectrum and will not affect the results of our 
analysis. 

Figure ^ shows the SEDs of the objects in our Class 
II sample, while figure displays those of our Class III 
objects. For all objects with measured extinction, their 
fluxes have been corrected for reddening using the ob- 
ject' s Ay from Ta bles |2 and [3] and Mathis's reddening 
law (|MathisHl99"o() with Ry = 3.1. Included are photo- 
metric data points from the optical to the mid-IR mostly 
from KH95, the 2MASS J, H, and K s fluxes, the IRAS 
12, 25, an d 60 itm fluxes, the IRAC 3.6, 4.5, 5.8, and 8.0 
/im fluxes ijHartmann et al.ll2005al) . where available, and 
the IRS spectrum. 

Also shown is the stellar photosphere for those objects 
with known spectral types; photospheric fluxes were de- 
rived from colors of main-sequence stars given in KH95 
and normalized at the J-band flux of each object with 
a spectral type later than G, and at the V-band flux of 
stars with an earlier spectral type (AB Aur, V892 Tau, 
RY Tau, and SU Aur). This normalization is based on 
the fact that for K and early M stars, most of the photo- 
spheric radiation is emitted around the J-band, but for 
earlier spectral types the peak emission shifts to shorter 
wavelengths. For spectral types earlier than G, the U- 
band would be better suited for normalization, but since 
U-band measurements are often poor and very sensi- 
tive to extinction, the V-band is a better choice. For 
F04101+3103, which is of Al spectral type, the J-band 
was used for normalization, since it might have an excess 
in the V-band. 

Th e IRAS data shown in the SED plots are mainly 
from iWeaver fc Jonesl l)1992j) . who used co-added IRAS 
data to determine fluxes even for faint objects (as many 
Class III objects are), which usually fell just below the 
sensitivity limit of IRAS. As can be seen in the com- 
parison with our IRS dat a, the fluxes for most Cl ass III 
objects, as determined bv lWeaver fc Jonesl l|1992l) . were 
often overestimated, resulting in mid-IR excesses for ob- 
jects whose emission is actually close to photospheric over 
the entire infrar e d ran ge. For objects not included in 
IWeaver fc Jonesl 1^)92), such as objects with identifiers 
from the IRAS Point Source or Faint Source Catalogues, 
we used the IRA S fluxes from these catalogs, mostly from 
the latter one. 

3.2.1. Class II SEDs 

All Class II objects, even the ones defined as WTTSs 
based on the equivalent width of their Ha emission line, 
show a clear mid-IR excess, indicative of the presence of 
circumstellar disks. The majority of our Class II objects 
are classical T Tauri stars, either single or in multiple 
systems; since they show a pronounced accretion signa- 
ture, indicated by their Ha emission or UV excess, it is 
expected that the circumstellar material in their inner 
disks generates a substantial infrared excess. 

There are only two single WTTSs in our Class II sam- 
ple: CoKu Tau/4 and IQ Tau. Their SEDs are quite 
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different from each other; IQ Tau has a substantial ex- 
cess from the near- to the mid-IR, indicating that at 
least the disk's inner regions (out to a few AU) are still 
prominent, while CoKu Tau/4 has no excess up to about 
8 ^m. This would indicate that WTTS disks are not 
all characterized by inner disk clearings, which is antici- 
pated given that they are no longer accreting. However, 
CoKu Tau/4 might be the only true single WTTS in our 
sample of Class II objects; it is a so-called transitional 
disk, where the inner disk has been cleared out, but a 
prominent outer disk is still present ijForrest et al.ll2004l 
iD'Alessio et al.ll2005D . Even though IQ Tau is classified 
as a WTTS based on its Ha equivalent width (7.8 A; 
iHerbie & Bellll988|) . it has a U-band e xcess, which pos- 
sibly indicates that is it still accreting l)Hartmann et al.l 
1998 derive d a mass accretion r ate of ~ 3 x 10 -8 M Q 
yr _1 , while [White fc Ghezll200l| only derived an upper 
limit of ~ 5 x 10~ 9 M yr _1 but noted its red K-L 
color). Therefore, WTTSs with prominent inner disks 
seem to be rare, confirming previous results of a short 
transition time from the CTTS s tate to a diskless WTTS 
state ffi-g-ISimon fc P^rSl. 

When comparing the SEDs of single and of multiple 
T Tauri stars, there are no obvious systematic differ- 
ences. However, since the IRS spectral range is sensitive 
to disk emission from a few tenths to a few AU from 
the star, we only expect to see any effects of multiplicity 
in close multiple systems (< 0'.'5, or 70 AU at the dis- 
tance of Taurus). Disks can surround each component 
or the entire multiple system; companions will truncate 
both circumstellar disks (at outer radii 0.2-0.5 times the 
semimajor axis of the binary) and circumbinary disks 
(at inner radii 2-3 times the semimajor axis of the bi- 
nary) , depending on the prop erties of the binary system 
ijArtvmowicz fc Lubowlll994j) . In wider systems, inner 
disks are not influenced by the presence of the other com- 
ponents and are thus expected to have the same proper- 
ties and lifetimes as disks around single stars. 

For most close multiple systems, the SED shows the 
combined photometry and spectrum of all components, 
and therefore the mid-IR flux could be generated by cir- 
cumstellar dust from disks around each component or 
from a disk surrounding the multiple system. Given that 
most binaries consist of CTTS pairs, both components 
will likely contribute to the IRS spectrum. 

By using ground-based near-IR photometry (K-band 
or, where av ailable, L-band) that resolved both com- 
ponents (e.g.. iWhite fc Ohe3l2001j) . we estimated each 
object's contribution to the IRS spectrum (see Table 
There is no clear correlation with binary separa- 
tion: while in some close (< 1") systems each compo- 
nent contributes roughly to equal amounts to the IRS 
spectrum (e.g., GN Tau), in some other close systems 
one component is the dominating one (e.g., V955 Tau). 
In the few mixed-pair binaries of our sample, where one 
component is a CTTS and the other a WTTS, the pri- 
mary component is usually the CTTS and dominates the 
mid-IR spectrum. Therefore, the IRS spectrum shows 
emission mostly from the circumprimary disk. This im- 
plies that the disk around the CTTS is either preferen- 
tially replenished or was formed as a more massive disk; 
the companion ends up with a smaller disk and a lower- 
mass star. This would also indicate that each system is 



strongly shaped by its initial conditions, i.e. the mass of 
each cloud fragment that formed a binary component, 
its kinematics, and whether other components formed 
nearby. 

Almost all of the Class II objects display pronounced 
10- and 20-/im emission features due to silicates. Usu- 
ally the 10-/xm feature is stronger than the 20- (im feature, 
and it also varies from object to object due to different 
grain composition and sizes. A smooth 10-/^m feature 
centered around 9.8 /jm is indicative of amorphous sili- 
cates (similar to the interstellar medium), while a struc- 
tured 10-/xm feature (see, e.g., the IRS spectrum of DH 
Tau) reveals the presence of crystalline silicates, likely 
processed fro m amorphous silicates in the circumst ellar 
environment l|Sargent et all200alFbrrest et alJ2004]) . In 
a separate paper, we explore the correlations between 
disk properties and the structure and strength of the 
10-/im silicate feature (D. M. Watson et al. 2006, in 
preparation) . 

In addition to the almost ubiquitous presence of sili- 
cate emission features, we note the close to complete ab- 
sence of PAH emission features in our objects of spectral 
type later than about Gl, with the possible exception of 
UX Tau A (see section l3~3*|) . PAH emission is typical for 
the more massive counterp arts of T Tauri sta rs, the Her- 
big Ae/Be stars (see, e.g.. lMeeus et al J 1200 1V and their 
absence in the circumstellar environment of lower mass 
stars is probably a result of the weaker UV radiation field 
of these stars. 

3.2.2. Class III SEDs 

All Class III objects in our sample have narrow Ha 
equivalent widths - all of them are WTTSs - and are 
not accreting within the detection limits, indicative for 
dispersal of their inner disks. Our IRS spectra support 
this finding: almost all Class III mid-IR spectra are pho- 
tosphcric. Even though the dispersal of disk material 
would suggest that these objects are more evolved, and 
thus older, the ages of Class III objects do not seem 
to be systematically larger than those of Class II ob- 
jects, which are supposedly in an ear lier phase of evo- 
lution (see, e.g.. I White fc Ghedl200l|) . It appears that 
the timescales for circumstellar disk evolution are rather 
determined by the initial conditions and the specific envi- 
ronment around a young star than set by some universal 
processes. Young Class III objects (< 1-2 Myr) appar- 
ently dispersed their circumstellar material faster than 
some long-lived Class II objects; m ultiplicity could play 
a role (see, e.g.. iMever et al.lll997|) but is likely not the 
only factor. 

Three of our Class III objects have nearby companions 
with infrared excesses that partly entered the 10'.'6 wide 
LL slit but not the narrower SL slit: DI Tau, HBC 423, 
and HP Tau/G2. DI Tau is separated by 15" from DH 
Tau, which is a classical T Tauri star (see figure^for the 
SED of DH Tau). In the case of HBC 423, V955 Tau, 
which is 10" away, and, to a small degree, HBC 422, 
which is 26" away, contribute to the flux past 14 fim (see 
figure [2 for the SED of V955 Tau). HP Tau/G2 is 10" 
from HP Tau/G3; both stars contribute to the LL flux, 
and in addition light of HP Tau (see figure ^) , which is 
21" from HP Tau/G2, probably entered in LL, especially 
toward the longer wavelengths. Given that the fluxes of 
DI Tau, HBC 423, and HP Tau/G2 are at photospheric 
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levels over the SL wavelength range (5-14 /im), but show 
a sharp increase where LL starts (at 14 /im), it is very 
likely that all of their apparent excess is caused by their 
companion(s). For these sources, we only show their SL 
spectra. 

Two more Class III objects show a similar behavior to 
the three previously mentioned objects: HBC 427 and 
V819 Tau. These two stars do not have any known com- 
panions with infrared excesses; however, for each of them 
the 2MASS K-band image shows a nearby source that 
likely contaminated the LL observation. V819 Tau is 
separated by about 10" from a faint source to the south, 
whose flux in the K-band is about 30 times smaller than 
that of V819 Tau. HBC 427 has a "companion" about 
15" to the southeast, which is about half as bright as 
HBC 427 in the K-band. Given that either "companion" 
did not fully enter the LL slit, but still likely contributed 
to the flux in this module due to the sharp rise in flux at 
the SL-LL boundary, they might have substantial mid-IR 
excesses. Since the mid-IR flux of those nearby sources 
is not known, we cannot determine which fraction of the 
infrared excess we observe is due to the target itself; the 
almost photospheric flux levels of HBC 427 and V819 
Tau at wavelengths shorter than 14 /im indicates that 
the intrinsic infrared excess is likely small. 

V410 X-ray 3 seems to have an excess at wavelengths 
longer than 20 /im, but since its spectrum is very faint in 
that wavelength region (~ 2 mJy), it could be an artifact 
of the extraction. This object is one of the lowest-mass 
T Tauri stars in our sample; i ts spectrum was first pre- 
sented in iFurlan et al.l l)2005aj) . 

HBC 356, Hubble 4, and L1551-51 were too faint 
beyond 14 /im, so only their SL spectra could be ex- 
tracted. For HBC 392, LkCa 1, and LkCa 21, the longer- 
wavelength part of LL (beyond 21 /im) was too faint, so 
their spectra include only a part of the LL spectrum. 

3.3. Notes on Individual Objects 

3.3.1. 04303+2240, 04370+2559, 04385+2550, 
CoKu Tau/3, FV Tau, Haro 6-13, MHO 3, 
V410 Anon 13 

These objects are all highly reddened T Tauri stars 
with an Ay larger than about 5 mag (see Table EJ; 
in addition, 04303+2240 a nd Haro 6-13 have relativel y 
high mass accretion rates l|White fc Hillenbrand! [2004). 
04370+2559 has a comp anion that is likely a substellar 
object (jltoh et alJll'9^al200l . 04385+2550, also known 
as Haro 6-33, has been considered as a Class I object 
by some authors bas ed on its SEP shape and bolomet- 
ric te mperature fe.g.. lMotte fc A ndre 2001; Youn get al.l 
2003). It is detected at sub-mm and mm wavelengths and 
is somewhat extende d in the former, but not in the lat- 
ter w avelength range (|Motte fc Andrel2001l:lYoung et all 
2003). On the other hand, the high extinction could in- 
dicate that 04385+2550 is a Class II object seen edge-on. 
V410 Anon 13 is a very low-mass classical T Tauri star, 
whose I RS spectrum was already studied in lFurlan et al.l 
(20053. 

The objects 04370+2559, CoKu Tau/3, and FV Tau 
have a CO2 ice absorption feature at 15.2 /im in their 
spectra. Since the lack of other ice features, the presence 
of a prominent silicate emission feature, and the SED 
shape suggest that these objects are T Tauri stars and 



not protostars surrounded by an envelope, the CO2 ice 
absorption feature is likely due to dense cold material 
along the line of sight to these objects. The relatively 
high extinction towards all of these objects would sup- 
port the idea of molecular cloud material causing the 
CO2 ice feature. 

3.3.2. CZ Tau B, DG Tau, DR Tau, FS Tau (Aa+Ab), 
HN Tau A, XZ Tau B 

According to iWhite fc Ohed" l|2001j) . these T Tauri 
stars have veiled optical spectra, unusually red K - L col- 
ors (> 1.4 mag), and, where measured, large Ha equiva- 
lent widths, indicative of high accretion rates. Therefore, 
their luminosities and temperatures are uncertain, and 
thus extinctions, ages, and mass accretion rates are not 
well determined. These stars might be experiencing an 
episodic phase of high accretion, which c ould be accom- 
panie d by a larger amount of extinction ijWhite fc Ghed 
2001). 

CZ Tau. CZ Tau has a peculiar IRS spectrum, which 
has a positive slope from 5 to 20 /im, but falls sharply 
beyo nd 20 /im. Based on their K- and L-band flux ra- 
tios l|White fc Ghed 120011) . the two components of this 
0'.'32 binary are probably contributing in equal measure 
to the IRS spectrum, but the B component is in a high- 
accretion state and therefore likely suffering from larger 
extinction. This could explain the shorter-wavelength 
slope of the IRS spectrum, but not the longer- wavelength 
drop; for the latter, disk truncation might be responsible. 

DG Tau, DR Tau. These two objects are the only 
single stars in this sample of high-accretion-rate stars; 
04303+2240 could be added to this subgroup, too. It 
seems that for these objects the excess emission from the 
accretion shock strongly veils the emission of the pho- 
tosphere, possibly dominating the emission from the ul- 
travi olet to the optical (see, e.g., [Gullbring et al J 11993. 
2000). Since it is difficult to determine the extinction and 
the flux level of the photosphere, we did not add stellar 
photospheres in their SED plots; we also did not attempt 
to adjust the extinction, but instead adopted the values 
published in the literature. 

DG Tau is a "flat-spectrum" T Tauri star, i.e. it has 
a flat SED over the infrared spectral range. Like other 
"flat-spectrum" sourc es, it is likely surro u nded by an in- 
falling envelope fe.g.. lCalvet et al]ll994ft . iWooden et alJ 
(2000) reported variations from emission to absorption 
behavior in the 10-/iin silicate feature, in addition to 
changes in dust mineralogy, on the timescales of months 
to years. Our IRS spectrum, taken in 2004 March, shows 
a very weak 10- /im silicate emission feature whose shape 
indicates that crystalline silicates might be present. We 
note that the flux level of the IRS spectrum is much lower 
than the IRAS fluxes at 12 and 25 /im. This is partly 
an aperture effect, since the IRAS beam included flux 
contributions from DG Tau, DG Tau B (a Class I object 
about half as bright as DG Tau) , and FV Tau (a Class II 
source, 0.25-0.5 times fainter than DG Tau), while the 
IRS observed only DG Tau (with slit sizes of 3"6, 4"7, 
and ll'.'l). Part of the discrepancy could also be caused 
by intrinsic variability of the target (s). 

FS Tau. FS Tau A and B, which are separated 
by 20", are both surrounded by reflecti on nebulosi- 
ties (|Krist et al.l(l99l IWhite fc Ghedl2l)0T|) : in our IRS 
observations, only the sub-arcsecond binary FS Tau 
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Aa+Ab was included. Given the smaller aperture of the 
IRS, we note that the IRS fluxes are considerably lower 
than the IRAS fluxes of FS Tau. 

HN Tau and XZ Tau. In the HN Tau and XZ Tau 
binary systems, both components are classical T Tauri 
stars, but the one component with the high mass accre- 
tion rate is clearly dominating the IRS spectrum (HN 
Tau A and XZ Tau B, respectively). 

3.3.3. AB Aur, F04101+3103, V892 Tau 

These three objects are Herbig Ae/Be stars and there- 
fore more massive and brighter than classical T Tauri 
stars. For AB Aur, we adopted the spectral type 
(AO) a nd extinction (Ay=0.25 mag) from iDeWarf et all 
(2003) . F04101+ 3 103 w as assigned an Al spectral type 
bv iKenvon et alJ <I1990|); in its SEP plot, we used t he 
BVRI photometry of llvengar fc Parthasarathvl l)1997|) . 

The spectral type of V892 Tau published in the litera- 
ture varies from B9, with an extinction o f Ay=8.85 mag, 
to A6 and an extinction of 4 mag (e.s.. IStrom fc Stroml 
119941 IKenvon fc Hartmannl Il995() . Here, we adopted 
the earlier spectral type and higher extinction (ad- 
justed somewhat to Ay=8.0), as it yields a more re- 
alistic SED. The high extinction is also supported by 
evidence of a nearly edge-on disk around V892 Tau 
ijHaas. Leinert. fc Richichl ll997). In addition, for V892 
Tau we used the B-, V-, a nd R-band measurem ents from 
the Herbig-Bell Catalog llHerbig fc Belllll988L h e reafte r 
HBC), and the R-I color from IStrom fc Stroml |l994), 
since KH95 onl y list a V-band magnit ude. The R and 
I photometry of IStrom fc Stroml l)1994|) is about 1 mag 
brighter than the HBC measurements, probably a sign 
of variability or an aperture effect. 

In the SED plots, the photospheres of AB Aur and 
V892 Tau were normalized at V; both objects already 
have an excess in the J-band. Even though the near-IR 
excess suggests the presence of an inner disk around V892 
Tau, it is classified as a WTTS. This could be explained 
by its early spectral type, for which the identification 
as WTTS based on its narrow Ha em ission line profile 
might not be valid ijCalvet et alJl2004)) . 

We note that V892 Tau is a 60 mas binary, and 
it also has a faint T Tauri star companion 4'.'1 away 
(Smit het alJ l2005). Our IRS observations included all 
three stars, but the subarcsecond binary likely dominates 
the flux of the system. 

3.3.4. GG Tau, UX Tau, UZ Tau, V773 Tau 

All of these four objects form quadruple systems, the 
highest-order multiple systems in our sample. 

GG Tau. GG Tau is a hierarchical quadruple sys- 
tem, consisting of two pairs of binaries: components Aa 
and Ab are separated by 0'.'25, Ba and Bb by 1'.'48, and 
components A and B by ~ 10". Our observations were 
pointed at the A binary; while in LL some flux of the 
B components also entered the slit, they are so much 
fainter than A that their contribution is negligible. In 
fact, the GG Tau B components are of lat e M spectral 
type, and Bb is most likely a brown dwarf l) White et alJ 
1999). All four stars exhibit accretion signatures in the 
form of strong Ha emission and are therefore surrounded 
by accretion disks; additiona lly, GG Tau A is surrounded 
by a circumbinary ring fe.g.. lMcCabe. Duchene. fc Ghed 
2002). 



UX Tau. UX Tau consists of component A, separated 
by 2'.'6 from C, and by 5'.'9 from B, which itself is a sub- 
arcsecond binary. Both A and C entered the SL slit, 
and in LL all four components contributed to the IRS 
spectrum. However, the A component is likely to domi- 
nate the mid-IR spectrum; it is the only classical T Tauri 
star of the system (the other components are classified 
as WTTSs), and it is brighter than components B and 
C by a f actor of 7 and 24, re spectively, in the L-band 
(3.6 urm lWhite fc Ghedl200l(l . To only contribute sub- 
stantially beyond 13 /im, but not at 10 /im and shorter 
wavelengths, the other components would have to have 
very red spectra. 

The peculiar IRS spectrum of UX Tau suggests the 
presence of some material in the inner disk and a sub- 
stantial outer disk. Since the other stars are a few arcsec- 
onds away, they are unlikely to influence the disk around 
the primary. UX Tau A has a very weak 10-/im silicate 
emission feature, but a prominent 20-/im feature, as well 
as some small peaks beyond about 20 /im possibly indi- 
cating the presence of crystalline silicates. Interestingly, 
UX Tau A shows a PAH feature at 11.3 /im (forsterite has 
a feature centered at slightly shorter wavelengths that is 
wider than the feature we observe); the 6.2 and 7.7 /im 
PAH features could be hidden by the strong continuum. 
Since UX Tau A's spectral type is K5, and we usually do 
not observe any PAH features in stars later than G, it 
would be the latest spectral type in our sample to show 
any PAH features. This latter fact might partly be due 
to the weak silicate emission of UX Tau, which allows us 
to easily identify the 11.3 /im PAH feature. 

UZ Tau. In the UZ Tau system, component A, a 
spectroscopic binary, is separated by 3'/5 from the ~ 
0'.'4 binary Ba-Bb. All four components entered the 
IRS slits. UZ Tau A is surrounded by an accreting 
circumbinary disk; even though the B components are 
also surrounded by accretion disks, their outer disks 
are likely truncated llJensen Koerner. fc Mathieulll996l 
Simon. Dutrev. fc Guilloteaul 120001 lHartigan fc Kenvonl 
2003|) . Therefore, the A component likely dominates the 
IRS spectrum. We note that in the SED plot the IRS 
fluxes are below the IRAC data points, which are the 
sum of the measured fluxes of UZ Tau e (i.e., A) and 
UZ Tau w (i.e., Ba+Bb); this could indicate that the 
IRS observations were somewhat mispointed and there- 
fore missing part of the flux of the brighter component 
(which is ~ 5-7 times brighter than the B components 
in the IRAC bands). 

V773 Tau. V773 Tau is a very tight quadruple sys- 
tem: all four stars are found within 0'.'25 from each other. 
V773 Tau AB is a spectroscopic binary and a WTTS, 
V773 Tau C is a CTTS, and the fourth component, dis- 
covered bv iDuchene et alJ l)2003|) . is an optically faint, 
and therefore deeply embedded, "infrared companion". 
V773 Tau D gradually increases in brightness from 2 /urn, 
where it is the faintest object in t he system, to 4.7 urn , 
where it is the brightest source Ipuchene et all H)03) . 
Thus, component D is likely the dominating source in 
the IRS spectrum. 

The SED of V773 Tau is somewhat similar to that of 
V807 Tau, which consists of a primary CTTS and a close 
binary WTTS, with separations of less than 0'.'5 between 
the components: both systems have little excess emission 
in the near-IR, a very weak 10-/mi silicate feature, but 
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a prominent 20-/im feature, suggesting that their inner 
disks have been partially cleared, possibly due to the 
action of the close components in these systems. 

3.3.5. HK Tau 

HK Tau consists of two components, separated by 2'.' 3; 
while the primary is surrounded by a face-on disk, the 
secon dary's disk is oriented edge-on l)Stapelfeldt et all 
1998). We were not able to resolve the binary in our 
IRS observations, but the primary is likely much brighter 
than the secondary ( flux ratios of 21 and 30 in the K- and 
L-band, respectively: iWhite fc G"he3l2001() . 

3.3.6. RY Tau and SU Aur 

These two stars are the only early G-type T Tauri stars 
in our sample. As with the SED plots of the Herbig 
Ac/Be stars, we normalized their photospheres at V. The 
IRS spectra of both objects display pronounced 10- and 
20-/zm silicate features; the smoothness of the 10- /an fea- 
ture indicates that the dust is mostly amorphous. 

SU Aur has an infrared excess already in the J-band, 
which would point to the presence of inner disk regions. 
Previous determinations of its EW[Ha| of 2 A i ndicated 
that it likely is a WTTS ijKenvon et al.lll998j) . in ap- 
parent conflict with the existence of an inner disk; how- 
ever, given its Gl spectral type, the classification as a 
CTTS or WTTS based on EWfH al might not be appro- 
pri ate (see IWhite fe Bas"rll2003|) . Recent measurements 
by iGalvet et alJ l)2004() suggest that it has a wide Ha 
equivalent line width, typical of CTTSs. Also, its UV 
excess reveals that it is accreting; it has an inferred mas s 
accretion rate of - 5xl0~ 9 M Q yr _1 ijCalvet et alJ2004|) . 
Therefore, SU Aur is probably a classical T Tauri star. 
However, it has a peculiarity in its spectrum: it displays 
PAH emission features like those seen in Herbig Ae/Be 
stars lIMeeus et abll200Tl lAcke fc van den Anc keH l2004t 
iSloan et al.ll2005l) but not in the late-type young stars of 
our sample. 

3.3.7. T Tau 

T Tau is a triple system; the northern and southern 
component are separated by 0'.'7, while the southern bi- 
nary has an orbital separation of Of.'l. The northern com- 
ponent, which is more massive and the source detected 
in the optical, is not very veiled, while the southern 
components are heavily extinguishe d and might there- 
fore be surrounded by an envelope l|Galvet et alJ 11991: 
iKoresko et al.l l 1997). In the infrared, the southern com- 
ponents vary in brightness due to variable extinction; 
in addition, the brightness ratio of the two components 
is changing on timescales of several months to a year 
(Bec ket alJl2004j) . 

T Tau N and S probably contribute about equally to 
the mid-infrared flux, even though the southern compo- 
nents could be brighter espec ially at the longer wave- 
lengths (see iGhez et all 1199 The IRS spectrum of 
the T Tau system is the only spectrum in our sample 
of T Tauri stars that shows a 10-/im absorption fea- 
ture. Since it was observed before that the northern 
component likely has a 10-/im emission feature, while the 
southern c omponent displays a silicate absorption fe ature 
at 10 /Ltm (|Ghez et al.lll99lUVan Cleve et al.lll994|) . the 
fact that we observe a 10-/xm absorption feature in the 



spectrum of the T Tau system suggests that T Tau S 
is the brighter source already at 10 /im and is indeed 
surrounded by an envelope. 

3.3.8. FP Tau, V836 Tau, VY Tau, ZZ Tau 

These four objects were classified as Class III objects in 
KH95. However, our IRS spectra clearly show that they 
are Class II objects; their SEDs show excess above pho- 
tospheric flux levels roughly starting at the L-band (3.5 
/im) - V836 Tau, whose L-band flux is likely erroneous, 
already has some excess at K (2.2 /im) - and their SED 
slopes are shallower than expected for Class III objects. 
The lack of near-IR excess emission in FP Tau, VY Tau, 
and ZZ Tau could be explained by cleared-out inner disk 
regions and/or very low mass accretion rates. In fact, 
the sp ectrum of ZZ Tau is rem iniscent of the spectrum of 
St 34 (jHartmann et alJfeOORbt) . where a spectroscopic bi- 
nary caused partial clearing of the inner disk regions, and 
dust growth and settling are likely responsible for weak 
silicate features and low dust continuum levels. ZZ Tau 
is also a close binary, with a separation of only 0'.'04 be- 
tween the two components. Like St 34, it has little excess 
emission below about 8 /im, but it does have a prominent 
silicate feature with a narro w peak at about 9.2 um., mo st 
likely due to silica (see, e.g.. lBowev fc Adamsoni r2002). 

The dust ex cess of VY Tau is particularly interesting. 
iHerbid l|1977fl emphasized that VY Tau has had a re- 
markable history of outbursts, varying in photographic 
magnitude from a lower base of m pg — 14 up to — 10 
or brighter, on timescales of a year or less. These inter- 
mittent, highly-irregular outbursts may then plausibly 
be attributed to events of accretion from a disk, possibly 
mostly cleared in its inner regions. This object should 
be monitored more extensively for further insights into 
its behavior. 

3.4. Morphological Sequence 

We arranged the IRS spectra of our 85 Class II objects 
into a morphological sequence according to the shape of 
their SED and the strength of their 10- and 20- /im silicate 
features (see Figures 13 to EJ). In each figure the spectra 
are ordered so that objects with stronger silicate features 
are at the top and objects with a steeper SED and weaker 
silicate feature are at the bottom; however, this sequence 
is not strictly followed to prevent overlaps of spectra in 
the figures. In addition, the spectra in figures |H1 and |H1 
should be considered as "outliers", whose peculiar IRS 
spectra cannot easily be assigned to any group of our 
sequence. To avoid erroneous reddening corrections due 
to uncertain extinctions Ay, which could result in more 
pronounced 10- /im features for large values of Ay, figures 
01 to El show the IRS spectra without any reddening cor- 
rection. Since few T Tauri stars in our sample have large 
extinctions, most spectra show the actual shape of the 
SED in this wavelength range even without reddening 
correction. 

Our first group of objects (Fig. [3J, group A, shows 
pronounced silicate features and a flat or somewhat de- 
creasing SED beyond about 20 fjm. The next group (Fig. 
0J), group B, is very similar to group A but has somewhat 
weaker silicate features. Group C (Fig. 03) is also charac- 
terized by more or less prominent silicate features, but it 
has clearly decreasing SEDs beyond 18 /im. Among these 
objects, DH Tau has a somewhat peculiar SED: the slope 
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of its spectrum between 5 and 8 /im and between 20 and 
35 /im matches that of the other spectra in this group, 
while there is a steep increase in flux between 14 and 
20 /Ltm. Also, its silicate feature is relatively wide and 
sho ws substructure cha racteristic of crystalline silicates 
(see iSargent et aTll2006ft . 

The spectra in groups A, B, and C display a large 
variety in their silicate feature, which does not seem to be 
correlated with the slope of the spectrum. Some 10-/im 
features are smooth and narrow, indicating amorphous 
silicates (e.g., UY Aur), while other silicate features are 
wider and show distinct peaks, which is a sign of grain 
growth and the presence of crystalline silicates (e.g., DK 
Tau). 

The fourth group (Fig. |BJ), group D, is characterized 
by a weak 10 /im feature and an overall negative SED 
slope over the IRS range. Many weak silicate features 
are relatively square and wide and show characteristic 
peaks of crystalline silicates; trends between the relative 
strengths of crystalline and amorphous silicates and the 
FWHM of the silicate feature will be explored by D. M. 
Watson et al. (2006, in preparation). In some spectra 
of group D, the 20- /jm silicate feature is somewhat more 
prominent than the feature at shorter wavelengths. This 
trend becomes more pronounced in our group E (Fig. |7|) , 
where almost no discernible 10-/mi feature is left in the 
spectrum (F04570+2520 being the exception), and the 
SED is decreasing steeply. 

Figures |H1 and show spectra that are not part of our 
morphological sequence. The spectra in figure |H1 dis- 
play an overall rising SED over the IRS spectral range 
and prominent 10-/im emission features. V892 Tau, 
F04101+3103, and SU Aur are of earlier spectral types 
than most of the Class II objects in our sample and, as 
opposed to T Tauri stars of spectral type later than Gl, 
have prominent PAH features in their spectra. CoKu 
Tau / 4 , DM Tau, and GM Aur are so-called transitional 
disks ilRice et al.ll2003t iForrest et al.ll2004t iQuillen et all 
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terized by a lack or decrease of excess emission short- 
ward of about 8 fj.rn due to inner disk clearings and a 
steep SED beyond 13 /im. UX Tau A would seem like 
a transitional disk due to its steep increase in flux at 
about 13 /im, but it has a substantial excess above pho- 
tosphcric flux levels beyond 2 /jm and only a very weak 
10-/im silicate emission feature. As mentioned in section 
13.31 UX Tau is a quadruple system, but the A component 
is likely the dominating source; some disk evolution has 
likely occurred in this system. 

Figure El shows a few sources with roughly flat SEDs 
from 5 to 35 /nm, as well as the spectrum of CZ 
Tau, whose SED is peculiar, unlike those of any other 
object in our sample (see section 13. 3|) . The flat- 
spectrum sources in this figure are objects associated 
with Hcrb ig-Haro jets and nebulosities in their surr ound- 
ings l)White fc Hillenbrandll2004UKrist et alJll998t and 
references therein). These objects are probably in tran- 
sition between Class I and II stage, when their envelopes 
have not fully dissipated yet. These "young" T Tauri 
stars include, for example, DG Tau, T Tau, and XZ Tau. 

4. ANALYSIS 
4.1. Dust Growth and Settling 



The mid-IR excess emission from Class II objects 
can be divided into optically thin emission from sub- 
micron sized dust grains in the flared disk surface layer 
(i.e., the upper layers of the disk atmosphere), whose 
most prominent signature is the silicate emission fea- 
ture at 10 /Ltm, and a continuum component, which orig- 
inates from deeper, optically thick layers of the disk at- 
mosphere. As disks evolve in ti me, dust grains grow 
and settle towards the midp lane l|Weidenschillingl|l997t 
IDullemond fc Dominikll2004|) : the upper disk layers close 
to the star (within a few AU) are depleted on a very 
short timescale, muc h less than the ~ 10 Myr lifetim e of a 
protoplanetary disk l|Dullemond fc Dominild l2"004) , with 
larger grains settling faster than the small grains. The 
larger grains around the midplane of the disk generate 
more emission at mm wavelengths, while the depletion 
of small grains in the upper disk layers causes a decrease 
in the continuum emission from the mid- to the far-IR, 
but these small grains still generate an emission band at 
10 /Ltm. Therefore, the effects of dust settling, and thus 
the extent of di sk evolution, can be seen clearly at mid- 
IR wavelengths i jMivake fc N akagaw al 19951 iFurlan et alJ 
l2005ttnTAlessio et a.l .1120061) . 

However, the idea of dust evolving from small grains 
to larger grains and eventually to planetesimals over the 
course of a few Myr is probab ly too simplistic: recently, 
IDullemond fc Dominild l)2005fl found that dust growth by 
coagulation occurs on very short timescales (less than 
~ 1 Myr, the typical age of a T Tauri star) in their 
simulations, and argued that small grains must therefore 
be replenished to produce the observed infrared excess. If 
grain growth and fragmentation both occur at the same 
time, grain sizes cannot be used to infer the evolutionary 
state of a disk, i.e. the presence of large grains in disks 
does not necessarily mean that a disk is older and thus 
more evolved. 

Over the last several years, it has become increasingly 
clear that disk models with dust settling are needed to 
explain the observed SEDs of T Tauri stars: maximum 
grain sizes and dust-to-gas mass ratios are different for 
the disk interior and the disk surface, with larger grains 
(< 1 mm) close to the midplane and small, ISM-like dust 
in the disk surface layer. To quantify the amount of 
dust settling, we use a variable, e, which is the ratio of 
the dust-to-gas mass ratio in the disk atmosphere and 
the standard du st-to-gas mass ratio o f the interstellar 
medium (1:100) l|D'Alessio et al.ll200l . A value of e=l 
means that no settling has taken place, while a small 
value of e indicates the presence of larger grains towards 
the midplane of the disk and fewer small grains in the 
upper layers. 

The morphological sequence from group A to E could 
be partly understood in terms of disk evolution: as dust 
gro ws, the 10- /xm silicate fe ature becomes wider and flat- 
ter ijPrzvgodda et alj|2003l) . and as it settles, the mid-IR 
excess decreases, causing a steeper negative slo pe of the 
mid-IR part of the SED ijD'Alessio et al.ll2006D . Due to 
small dust grains still present in the disk atmosphere, 
the 10-fim feature can still be strong relative to the con- 
tinuum; it is more pronounced than the 20-/tm feature 
due to the wavelength-dependent emissivity of the dust. 
A decrease in the 10- to 20-/im flux ratio could indi- 
cate a decrease of emission from inner disk regions rela- 
tive to the outer parts, since the 10-/Ltm emission arises 
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mostly from inner regions of the disk, while the main 
20-/im emission region lies further out. Given that dust 
growt h and settling is fas t est in the inner parts of the 
disk ijHavashi et al.l 119851: iWeidenschillind 1 1 997(1 . it is 
expected that the silicate emission from the inner disk 
would decrease first. A more detailed, quantitative study, 
which is beyond the scope of this paper, is required to 
tie our morphological sequence to disk evolution; here, 
we tentatively suggest that dust has grown and settled 
to a larger extent in the disks of groups D and E than in 
groups A, B, and C. 

4.2. Settled Disk Models 

To analyze our Class II objects in a more quantitative 
way, we compare them to a grid of accretion disk models 
that include dust settling l(D'Alessio et alJ 12006). The 
vertical disk structure is de r ived s elf-consistently, as in 
iD'Alessio et"aTI IpM Il999l l200lh : the resulting flared 
disk extends from the dust sublimation radius to an outer 
radius of 300 AU, which was set arbitrarily, but does not 
affect the shape of the mid-infrared spectrum. At the 
inner disk edge, the "wall" , light from the star is incident 
parallel to the surface normal. The disk is heated by 
viscous dissipation and stellar irradiation. The central 
source is a typical T Tauri star in Taurus: T=4070 K, 
R=1.864 Rq, M=0.8 M q , L=0.79 L q . 

Silicates and graphite mak e up the dust grains in the 
disk (optical constants from iDraine & LeejTl98l : their 
sizes range from 0.005 to 0.25 /im in the upper disk layers 
(i.e., similar to the interstellar medium), while in the 
lower disk layers, the maximum grain size is set at 1 mm. 
The dust particles follow a size distribution of the form 
n(a)da — a~ 3 5 da between the minimum and maximum 
grain sizes. Models with either isotropic scattering of 
the incident radiation by dust grains or perfect forward 
scattering were considered, with the true case lying in 
between the two limiting cases. The value of e was also 
varied; to include the effect of both dust settling and 
grain growth, a higher depletion of small grains in the 
upper layers is accompanied by an increase of large grains 
in the lower disk layers. Thus, while the dust is depleted 
in the disk atmosphere, the amount of dust is increased 
in the layers around the midplane, conserving the total 
amount of dust. The models of the grid vary in accretion 
rate (10 -8 and 10 -9 M Q yr _1 ), inclination angle along 
the line of sight (75?5, 60°, 50°, 40°, 30°, 20°, and 11?5), 
amount of dust settling (e— 1, 0.1 0.01, and 0.001), an d 
dust composition. As we did in iFurlan et all l|2005b|) . 
we use the spectral coverage of the IRS combined with 
models to estimate the degree of dust settling. 

4.3. Spectral Indices 

4.3.1. Definition 

To characterize the SEDs of the T Tauri stars in our 
sample, we computed spectral indices at different wave- 
lengths in the IRS spectral range that r epresent con- 
tinuum emission, similar to our analysis in lFurlan et al.l 
( 2005bj). We integrated the flux of our dereddened spec- 
tra in the following wave bands: 5.4-6.0 /im (center wave- 
length, A c , at 5.7 /im), 12.5-14.0 fim (A c = 13.25 /im), 
and 23.5-26.5 /im (A c = 25 /im). Then we divided the 
result by the width of the respective wave band, which 
yielded a flux density in wavelength units (F\) for each 



band. To obtain spectral indices, n, we computed 



log 



A^J /log u 



first between 5.7 and 13.25 /im (^6-13)1 then between 
13.25 and 25 /im (^.13-25)- 

Furthermore, to quantify the properties of our morpho- 
logical sequence, we used our IRS spectra not corrected 
for reddening and calculated a spectral index covering 
most of the IRS spectral range, from 6 to 25 /im (77,6-25), 
and estimated the strength of the 10-/im silicate feature. 
For the latter, we interpolated the continuum from data 
between 5.0 and 7.5 /im and between 13.0 and 16.0 /im 
with a third-order polynomial, subtracted it from the sil- 
icate feature defined between 8.0 and 12.4 /im, integrated 
the flux of the continuum-subtracted silicate feature, and 
normalized it to the continuum integrated over the same 
wavelength range. Due to uncertainties in the contin- 
uum fit, this procedure might not always yield accurate 
measurements of the 10- /im silicate feature strength, but 
for the majority of objects it will result in reasonable 
estimates. 

4.3.2. Results: IRS Data 

Figure ^| displays the continuum-subtracted, inte- 
grated flux of the 10-/tm feature, normalized to the con- 
tinuum, versus the 71.6-25 spectral index. The data points 
belonging to spectra from groups A to E of our mor- 
phological sequence are identified by different plotting 
symbols; the data for the plot can be found in Table 
0] As indicated by our morphological sequence, ob- 
jects with the smallest values of 71,6-25 in our sample 
(n 6 _25 < —1.0) have weak 10-/im silicate features. The 
spectra of groups D and E have the smallest 10-/im fea- 
ture strengths, but SED slopes in the range of those of 
group C, which has stronger silicate features. Only ob- 
jects with ng-25 > —0.5 have very strong features, but 
for these values of the spectral index there is a large 
spread in silicate feature strength. 

Figure ITT1 shows a plot of 77,13-25 versus tt-6-13, com- 
puted using the IRS spectra of the Class II and III ob- 
jects in our sample (open diamonds and crossed squares, 
respectively), and for two extremes, a flat, passive disk 
and a naked photosphere; the data for the plot are also 
given in Table A photosphere in the mid-IR has a 
spectral index of —3 (Rayleigh- Jeans limit), while a geo- 
metrically thin, optically thick disk has an index of —4/3. 
Due to some spectrophotometric uncertainty in our IRS 
data, the typical error bars for the spectral indices are 
about ±0.05. The dashed boxes delineate outliers and 
will be discus sed below. 

As noted in lFurlan et all l|2005b|) . the fact the Class II 
objects in Fig. ^]have spectral indices that are larger 
than the value expected for a flat, optically thick disk 
means that their disks must be flared. Only one object 
in our sample, F04570+2520, has an 77,13-25 spectral in- 
dex that is slightly smaller than that of a flat, optically 
thick disk; its SED shows only a weak 10-/im silicate 
feature and a steep decrease in flux beyond it, probably 
indicating an evolved disk in which a large fraction of 
the dust has grown and settled. For most objects, the 
n6_i3 and ?ii3_25 spectral indices seem to be correlated; 
a steeper slope of the spectrum between 6 and 13 /im is 
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usually accompanied by a steeper slope from 13 to 25 /im. 
As can be seen from the frequency distribution of both 
indices in Figure [T51 the distribution peaks at values be- 
tween — 1.5 and for ne-13, and between —1 and 0.5 for 
"■13-25- While negative spectral slopes larger than ~ — 1 
prevail for both wavelength ranges, which is typical for 
Class II objects, the slope measured between 13 and 25 
/im tends to be natter than the one measured between 6 
and 13 /im. 

In addition, there is a small number of objects with 
spectral indices 7113-25 larger than about 0.5, which is 
represented by the seemingly randomly scattered points 
at positive 7113-25 values in Figure ^2 an d framed by 
dotted boxes in the figure. These "outliers" can be di- 
vided up into three groupings: the 6 objects in the lower 
left box, the 4 objects in the right box, and the 4 objects 
in the uppermost box. 

The first-mentioned grouping contains DH Tau, LkCa 
15, AB Aur, HK Tau, Haro 6-13, and SU Aur. These 
Class II objects are characterized by a steep rise in slope 
starting at 13 /tm, indicating a strong 20-/tm silicate fea- 
ture, and a relatively flat SED past 20 /im. The right 
box frames F04101+3103, MHO 3, V892 Tau, and, at a 
negative 7113-25 value, CZ Tau. These objects also have a 
strong 20-/im silicate feature, in addition to a pronounced 
10-/tm emission feature, but their SEDs decrease past 20 
/im. Also, their emission from 5 to 8 /im is lower than 
the emission at longer wavelengths, resulting in a posi- 
tive index "6-13. The decrease of the SED beyond about 
20 /xm is especially dramatic for CZ Tau, which explains 
the low value of 7113-25 of this object. As mentioned ear- 
lier, CZ Tau's outer disk might be truncated, resulting 
in a substantial decrease of its far-IR flux. 

The final grouping includes the known transitional 
disks in Taurus, CoKu Tau/4, GM Aur, and DM Tau, as 
well as UX Tau A at "6-13 ~ —2. These "outliers" have 
negative spectral indices from 6 to 13 /im and positive 
spectral indices from 13 to 25 /mi, the latter due to a 
steep increase in flux at about 13 /im. 

Also included in Figure El ar e the spectral indices 
"13-25 versus «6-i3 f° r those Class III objects in our 
sample whose 25 /tm flux was detected. However, in 
many cases the 25-/im flux is very weak and the LL spec- 
trum noisy, resulting in a larger uncertainty for the index 
calculation. Nevertheless, most of the indices of our Class 
III objects are close to values of —3, as expected for close- 
to-photosphcric fluxes. One object, V410 Tau, which was 
not included in the plot, lies at spectral indices smaller 
than —3, which is mainly a result of its very noisy LL 
spectrum; its IRS spectrum seems to follow the slope of 
a photosphere, with a possible slight excess. V410 X-ray 
3 lies at 7113-25 = —1.6 and "6-13 = —2.9; HBC 427 and 
V819 Tau, which have an apparent infrared excess in LL, 
can be found at 7113-25 = —0.7 and n 6 -i3 = —2.8 and 
—2.6, respectively. If at least part of their infrared ex- 
cess is real, the latter three objects would constitute the 
transition phase between Class II and Class III, which 
seems to be very rapid due to the lack of objects with 
spectral indices between —1 and —2. This result is sim- 
ilar to what is found when the K -L colors are plotted 
versu s the K-N colors (see, e.g.. iKenvon fc Hartmaniil 
1995): there is a continuous transition between Class I 
and II objects, with the latter displaying bluer IR colors, 
but a gap between the Class II sources with the small- 



est infrared excess and the Class III sources, which have 
almost photospheric colors. 

4.3.3. Results: IRS Data and Models 

We computed analogous spectral indices (7113-25 vs. 
7ig_i3) for the grid of accretion disk models and over- 
plotted the results from IRS data and disk models in 
Figure El We show indices computed from two sets of 
disk models, one with isotropic scattering of the incident 
stellar radiation by dust grains and the other with per- 
fect forward scattering (which has the same effect as no 
scattering, i.e. albedo=0); they represent the two limit- 
ing cases of scattering. There is only a slight difference 
between these two sets; spectral indices are not very sen- 
sitive to the degree of scattering in the models. Each 
model set contains models computed with two mass ac- 
cretion rates, 10 -8 and 10~ 9 M Q yr -1 , which are typi- 
cal for T Tauri stars in Taur us (see, e.g. iGullbring et all 
19981 IHartmann et al.l lT998h. 

The spread in the 7i6-i3 index is partly due to differ- 
ent accretion rates, but mostly an effect of varying incli- 
nation angles. This can be understood in terms of the 
wall emission, which dominates the flux at wavelengths 
shortward of about 6 /im; the wall emission increases as 
the inclination angle increases, resulting in a larger in- 
crease of the 6-/im flux than of the 13-/im flux, while 
the disk emission, wh ich dominates at 25 /im, decreases 
Ip'Alessio et al.ll2006D . The different extents of dust set- 
tling in the models, which are described by e, cause a 
spread in both 716-13 an d 7113-25 . Since for mass accre- 
tion rates of 10~ 8 and 10~ 9 M yr _1 i rradiation heat- 
ing represe nts the main heating source l|D'Alessio et alJ 
1999, 2006), more settling causes a decrease in the mid- 
IR continuum emission; the disk intercepts less radiation 
from the central star if its outer parts are less flared. 
This affects the 25- /tm flux more than the flux at shorter 
wavelengths, resulting in smaller spectral indices (i.e., 
steeper SED slopes) for smaller e. 

For example, AA Tau, w hich has a mass accr etion 
rate of 6.5x 10" 9 M yr" 1 l|White fc Ohedl27)0l . lies 
at 7i6-i3 = —0.87 and 7113-25 = —0.37, which is close 
to the model with an accretion rate of 1 0~ 8 yr" 1 , 
e = .001, and inclination angle of 50°. iBouvier et alJ 
( 1999) derived a lower limit for the inclination angle of 
AA Tau of 53°, but favored i — 75° to explain the pe- 
culiar variability of this object. Our comparison would 
suggest a somewhat lower inclination angle, and in addi- 
tion a large amount of settling in the disk of AA Tau. As 
another ex ample, BP Tau, which has an inclination an- 
gle of 39° l)Muzerolle et al.l |20fl3 l bj)_ and a m ass accretion 
rate of 1.3x 10" 8 M vr" 1 l|White fc Chezll200l lies at 
7ig„i3 = —0.67 and 7113-25 = —0.24, which is in between 
the model points for M = 10~ 8 M yr^ 1 , e = 0.001, and 
inclination angles of 30° and 40°, respectively. 

Most Class II objects can be described with models 
that include substantial settling, with an e between 0.01 
and 0.001, which means dust depletions of factors of 100 
to 1000 in the upper layers. The spectral indices for ac- 
cretion disks with e < 1 are lower than for those in which 
no settling has taken place, since the reduced mid-IR con- 
tinuum causes a steeper slope, towards the value of —4/3 
of a flat, optically thick disk. The settling of dust does 
not affect the silicate emission features, which we see in 
virtually all of our Class II spectra; these features are 
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generated by the small grains still present in the upper 
hot layers of the disk atmosphe re. As indicated ea rlier 
with a smaller sample of CTTSs llFurlan et al.l2005bj) . we 
conclude with this larger sample of CTTSs and WTTSs 
that the majority of T Tauri stars have disks in which 
dust has started to grow and settle towards the mid- 
plane, thus affecting the SEDs of these young (1-2 Myr) 
pre-main-s equence stars. This is consisten t with the sim- 
ulations of Dullcmond & Dominik (2005), which predict 
dust growth on timescales of ~ 1 million years or less. 

Given that all of our T Tauri stars are at about the 
same age, the difference in the spectra of various objects 
must be due to different initial conditions and different 
circumstellar environments influenced, for example, by 
the presence of close companions. Thus, the extent to 
dust settling in a disk, at least in its initial stages, might 
be less a function of age and more related to the object's 
formation history and current environment. 

We caution that the comparison of our data with a 
grid of accretion disk models is valid to gauge trends in 
our data, but that additional parameters, in particular 
the composition of the dust, need adjustments before any 
detailed interpretations can be carried out reliably. For 
example, choosing glassy pyroxene or olivine instead of 
"astronomical" silicates and graphite will generate dif- 
ferent dust emission, both in the continuum and in the 
emission bands, and different spectral indices will result. 
Therefore, the detailed composition of the dust adds ad- 
ditional scatter to the spectral indices plot, which has to 
be considered in addition to the effects of inclination and 
dust settling. The values we derived for e should thus be 
treated as rough estimates only. The effect of dust com- 
position on the models, as well as models of individual 
objects, will be presented in a future paper (P. D'Alessio 
et al. 2006, in preparation). 

4.4. Taurus Median 

In order to determine a typical SED range for T Tauri 
stars in Taurus, we computed a median SED from 1.25 
to 34 (im using the 2MASS J, H, and K s photometry 
and the IRS spectra of all our 85 Class II objects, dered- 
dened using the Ay values listed in Tabled To obtain a 
clearer representation of the flux values in the IRS spec- 
tral range, we integrated the IRS spectra over narrow 
wave bands and divided by the width of the band, creat- 
ing the equivalent of narrowband photometric measure- 
ments at 5.7, 7.1, 8.0, 9.2, 9.8, 11.3, 12.3, 13.25, 16.25, 
18.0, 21.0, 25.0, 30.0, and 34.0 /mi. The width of the 
bands ranged from 0.6 /xm at the shorter wavelengths to 
3 fim at the longer wavelengths; the width of the sili- 
cate bands, cent ered at 9.8 and 18 /im , was chosen to 
be 2 /im. As in iD'Alessio et alJ (1999), we normalized 
all flux values at H (1.65 /im) before computing median 
values, since, for most objects, the H-band flux is pho- 
tospheric, and thus normalization at H is equivalent to 
normalizing to the stellar luminosity. Since for most of 
the T Tauri stars in our sample mass accretion rates are 
low enough (< 10~ 7 M Q yr _1 ) that the star is the main 
heating source of the disk, the infrared excess of the disk 
should scale with stellar luminosity. By computing the 
median values, we also obtained the quartiles, which de- 
fine the range around the median where 50% of all flux 
values lie. 

Figure 1141 shows our median with the quartiles (in- 



dicated by th e erro r bars) and that obtained by 
ID'Alessio et, all l)1999ft . We computed the median using 
all 85 Class II objects in our sample (upper, green data 
points), and by including only stars with spectral types 
between K5 and M2 (lower, red data points). The latter 
data set constitutes a more homogeneous sample of stars 
with similar effective temperatures and allows a more 
direct comparison to the median from ID'Alessio et alJ 
(1999); it is also given in table El 

With the IRS data, we are able to resolve the silicate 
emission features at 10 and 20 /zm; the median SED 
clearly shows that most objects have pronounced silicate 
emission and thus small grains in the hot upper layers 
of the disk atmospheres. It is remarkable that our me- 
dian SED at 25 iim a lmost coincides with the one from 
ID'Alessio et all ljl999ft . who used IRAS 25 /jm flux mea- 
surements. The median SED for K5-M2 stars is lower 
than the one for all Class II objects in our sample, as 
expected for a sample of stars that excludes earlier, and 
thus hotter, spectral types. In particular, the decrease 
in the 13-25 /im range is more pronounced than that at 
shorter wavelengths, indicating that the silicate feature 
at 20 /im decreases more than the 10-/mi silicate feature 
for late spectral types. 

Figure ^] displays the IRS median SED for spectral 
types K5-M2 (red), the med ian computed from IRAC 
data l|Hartmann et alJ l2005ai purple), and the photo- 
spheric component ( dotted line) , represented by th e SED 
of the WTTS HBC 427 ijHartmann et alJ l2005a). We 
compare these to SEDs of two irradiated accretion disk 
models with different degrees of dust settling. The pho- 
tospheric component is normalized at J, while all the 
other SEDs are normalized at H. 

The models shown in figure E3 we re constru cted ac- 
cording to the methods of D'Alessio et af] 1)2006(1 and are 
described in section 14.21 with a few differences: accre- 
tion shock irradiation is included as heating mechanism 
for the disk (but it does not make a significant differ- 
ence for M < 10~ 8 M Q yr _1 ), and the dust mixture con- 
sists of spherical particles of amorphous pyroxene (glassy 
Mgn 7 Fen s Si03, optical constants from lDors chner %\ alJ 
1995) and graphite (optical constants fromlDralne fc Lee} 
1984). This dust composition results in higher silicate 
emission at 10 and 20 /im and therefore is a better fit for 
the Taurus median SED (see also discussion in § I4.3.3fl . 
The models were computed for an inclination angle of 
60°, a mass accretion rate of 3 x 10~ 8 M Q yr _1 , and 
two values of the settling parameter e, 0.1 and 0.001. 
Comparison of the models with the median and quartiles 
suggests large degrees of depletion in the upper layers of 
disks in Taurus: about 50% of the disks have a dust de- 
pletion larger than a factor of 1000 (e < 0.001), while 
about 25% of the disks can be characterized with de- 
grees of depletion between a factor of 10 and 1000. Only 
about 25% of the disks have between a factor of 10 and 
no depletion of dust in the upper disk layers. 

4.5. Class III Objects: Upper Limits for Disk Masses 

We estimated upper limits for the mass of small, warm 
dust grains (size < 10 /im, temperature > 100 K) in the 
disks around the 26 Class III objects we observed, assum- 
ing optically thin dust emissi on from "astronomic al" sili- 
cates (optical constants from lDraine fc Leall98 4"). Since 
none of the Class III objects considered here show sili- 
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cate emission bands at 10 or 20 /im, larger grains (> 3-4 
/im), if any, must be present in their disks; for this rea- 
son we included only grains between sizes of 3 and 10 /im 
(the mid-infrared is not sensitive to emission from larger 
particles). First, we calculated the equilibrium tempera- 
ture for each grain size at a distance r from the star by 
solving the following equation numerically: 

i-f] I Qabs{y)B v {T*)dv = A J Q abs {v)B„{T d )dv, 

where i?„ is the stellar radius, Q a bs is the absorption 
efficiency, and B„(T) is the Planck function, evaluated 
at the stellar effective temperature T* and at the dust 
grain temperature T d , respectively. Assuming that the 
disk is axially symmetric and geometrically thin, we 
parametrized the optical depth as 

T~(r,a) = £ (r ) ( — ) (— ) Tra 2 Q abs (a, v), (2) 

where £o( r o) is the reference surface number density of 
dust grains of size ao at a distance ro from the star; this 
column density varies with radius as r -7 , and the grain 
size distribution varies as a -3 5 . We chose Tq = 1 AU, 
ao = 1 ^m, and, besides for a few exceptions, 7 = 1. The 
flux emitted from the disk at each wavelength is obtained 
by integrating over annuli and summing over each dust 
grain size: 

Fx = E (^2 l ^ 2nrT ^ a 05„(T d (r, <H))drj , (3) 

where r(r, Oj) is defined in equation P|l. 2d(r, flj) is ob- 
tained from equation lfT|l. and D is the distance to the 
star, assumed to be 140 pc, the distance to the Taurus 
star-forming region. The total dust mass in the disk is 
given by 

rTout 

M d = ^Tnid / 2irrE(r)dr, (4) 

where m,d is the mass of a single dust grain. To obtain 
upper limits for the dust masses in our Class III objects, 
we adjusted the parameter £o( r o) (and in some cases also 
7) , until the disk emission from equation J3J matched the 
emission of each Class III object after photosphere sub- 
traction. The photospheres are re presented by AMES- 
dusty models of lAllard et al.1 lj2001() for effective temper- 
atures up to 5000 K; for hotter stars, we approximated 
the photospheres by blackbodies in the Rayleigh- Jeans 
limit. We integrated the flux from 0.1 to 50 AU and 
summed over 11 grain sizes from 3 to 10 /im radius. 

With all the assumptions mentioned above, our model 
fits result in typical upper limits for warm, small dust 
grains in any circumstellar disks around our Class III ob- 
jects of a few times 10 lunar masses. Assuming that 
the infrared excess of HBC 427 and V819 Tau can be 
fully attributed to disks around these objects, the maxi- 
mum amount of dust mass in these disks would amount 
to a few times 10 -3 lunar masses. In addition, since their 
excess starts rather abruptly at 14 /im, their disks would 
require inner holes of a few AU in size to explain the 
lack of shorter- wavelength excess emission. HP Tau/G2 
seems to already have a slight excess shortward of 14 /im; 



again assuming all the excess emission originates in a 
disk around this object, the upper limit to the dust mass 
would amount to somewhat less than 10~ 4 lunar masses. 
A dust mass of a few times 10 -4 lunar masses is com- 
parab le to that found around debris disks ijHabing et alJ 
1200 1|) . although closer to the lower mass values. 

5. DISCUSSION 

The multi-wavelength SEDs, which weight the spec- 
trum in terms of the relative power emitted at each 
wavelength, provide an immediate initial estimate of 
the relative luminosities of the star and the disk. In 
most cases, the luminosity of the disk is less than that 
of the star because most of the disk heating and thus 
disk emission are due to absorption and reprocessing 
of stellar radiation. However, of the 85 Class II sys- 
tems shown in Figure [T] XZ Tau, UY Aur, HN Tau, 
FS Tau, 04187+1927, 04200+2759, DG Tau, and DP 
Tau appear to have anomalously large IRS excesses rel- 
ative to the central star. As noted in section 13.31 a 
few of these objects have envelopes; UY Aur has an 
IR companion separated by Of.' 9 from the primary and 
likely seen throug h the disk of the primar y, given its 
larger reddening IjHartigan &: Kenvonl 12003) . Together 
with the other ~ 9 objects with very high and there- 
fore uncertain Ay mentioned in section T3. 31 there might 
be about 17 candidates for edge-on disks. In addition, 
HK Tau B llStaoelfeldt et alJll998l) and like ly also AA 
Tau ijMenard et alJl2003UBouvier et alJl2003|) have edge- 
on disks, bringing the total of T Tauri stars likely seen 
at an inclination angle close to 90° to 19. Thus, 22% 
of the Class II objects in our sample (or 17% of the T 
Tauri stars of our Taurus sample) are likely seen edge- 
on, which is roughly what is expected from statistical 
considerations for objects with inclinations between 77° 
and 90°. In addition, these objects do not show any 
dee p silicate absorption fe atures, as expected from mod- 
els (jD'Alessio et al.lll999|) : therefore, either these highly 
inclined disks are not massive enough, or we are ob- 
serving the 10-/im thermal emission from the inner disk, 
scattered by the dust in the oute r disk atmosphere (see 
iMcCabe. Duchene. fc Ghe3l2003l) . 

A few of our Class II objects show a decrease or lack of 
IR excess at shorter IRS wavelengths, indicative of inner 
disk clearings or gaps. Some of these sources (GM Aur, 
DM Tau, CoKu Tau/4) have been descri bed as transition 
disks and are the subject of other papers l)D'Alessio et alJ 
l200fil ICalvet et al.ll2005h . Only 4% of the Class II ob- 
jects in our sample, or 3% of the HIT Tauri stars we 
observed in Taurus, are in this transitional disk phase, 
which would indicate that this stage is short-lived. 

If their longer-wavelength excess proves to be real, 
HBC 427, V819 Tau, and possibly also V410 X-ray 3, 
could also be considered as objects in transition, even 
though their excess emission is so weak that it is proba- 
bly not optically thick disk emission. On the other hand, 
they might be regarded as debris disks, where dust is gen- 
erated by collisions of larger bodies and is optically thin 
at all wavelengths. If that were the case, they would be 
among the youngest debris disks discovered so far. 

UX Tau A has a very peculiar spectrum, with signifi- 
cant excess at shorter wavelengths, but a sharp transition 
at 20 /im. Even though it is a multiple system, the other 
components are at distances that are unlikely to affect 
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the inner disk we are observing with the IRS. Also, their 
flux contribution in the mid-IR is expected to be small. 
Thus, the IRS spectrum of UX Tau A suggests that it 
might be in a pre-transition state, when most dust in the 
inner disk has already settled and grown, and thus the 
lifetime of the inner disk is coming to an end. 

Other objects that could be in a pre-transition state 
are 04385+2550, Haro 6-13, HK Tau, and LkCa 15; their 
20-/xm silicate feature is stronger or comparable to the 
10- /zm feature, reminiscent of transitional disks, but they 
have substantial excess emission in the near-IR. In these 
objects, a gap might be opening, but a significant inner 
disk is still present. 

6. CONCLUSIONS 

We presented the mid-infrared spectra of 85 Class II 
and 26 Class III objects in Taurus, obtained with the 
Infrared Spectrograph on board the Spitzer Space Tele- 
scope. This large sample of pre-main-sequence stars of 
similar age, belonging to the same star-forming region, 
allows us to draw the following conclusions about cir- 
cumstellar disk evolution: 

• All CTTSs and some WTTSs display mid-IR ex- 
cesses, which are generated in the inner (a few tenths to 
a few AU) parts of circumstellar disks. Accretion signa- 
tures like the Ha equivalent line width do not necessarily 
indicate the presence or absence of disks, but they usu- 
ally reveal whether the innermost disk regions are still 
present. While accreting T Tauri stars have dusty disks 
extending from the dust sublimation radius outward and 
thus generate excess emission starting at near-IR wave- 
lengths, WTTSs often have inner disk holes and thus 
lack near-IR excess emission. Since there are few WTTSs 
with disks, the transition period between an active (i.e., 
accreting) and a passive (i.e., purely reprocessing) disk 
stage seems short-lived. 

• The details of the mid-IR excess, and thus the distri- 
bution and properties of the dust in circumstellar disks, 
do not seem to be correlated with the multiplicity of the 
young star system. A sub-arcsecond (< 140 AU at the 
distance of Taurus) companion docs not necessarily in- 
fluence the inner disk; even spectroscopic binaries, like 
DQ Tau, seem to have substantial inner disks. We note 
that our IRS spectra cannot separate the contribution of 
close multiple systems, and therefore, in some cases, they 
show the combined mid-IR spectrum of circumprimary, 
circumsecondary, and circumbinary disk. The mid-IR 
spectra of multiple systems do not have any peculiarities 
that would set them apart from single stars; this could 
be explained if most disks in these systems are not in- 
fluenced by the presence of companions and either the 
circumprimary disk clearly dominates the system or all 
components have similar disks, or if a circumbinary disk 
mimicked the behavior of a circumstellar disk around a 
single star. 

• Almost all Class II objects show 10- and 20-ttm sili- 
cate emission features in their mid-IR spectra, generated 
in the optically thin disk surface layer. There is a wide 
variety of feature strengths and crystallinity, indepen- 
dent of the spectral type and multiplicity of the objects. 
However, objects with weak 10-ttm silicate features usu- 
ally have a decreasing SED in the mid-IR range, which 
could be an indication of dust growth and settling. In 
objects where the 10- and 20-ttm features are suppressed, 



dust grains must exceed radii of about 4 /zm. 

• After comparing the spectral indices 
d\og(\F\)/d\og(\) between 6 and 13 ttm and be- 
tween 13 and 25 ^m of our Class II objects with 
accretion disk models that include dust settling, we infer 
that our data are consistent with these models. Within 
the assumptions of the adopted models, most T Tauri 
stars seem to have experienced substantial dust settling 
(and likely also dust growth), with depletions of factors 
of about 100-1000 of the standard dust-to-gas mass ratio 
in the disk atmosphere. Thus, we are able to estimate 
the amount of dust settling for a large sample of T Tauri 
stars, using the high sensitivity and spectral coverage of 
the IRS. This result will aid in the refinement of models 
predicting the evolution of dust in protoplanetary disks. 

• Most WTTSs in our sample are Class III objects, 
which have nearly photospheric fluxes in the infrared, 
indicating a dispersal of dust in their circumstellar disks. 
The upper limit for the mass of small, warm dust grains 
in these disks lies around 10 -4 lunar masses, which is 
similar to the dust mass of low-mass debris disks. In 
some cases nearby companions (~ 10"), which seem to 
have an infrared excess, contaminated the LL part of 
our Class III sources. It is interesting that objects of 
roughly the same age, formed in close proximity to each 
other, do evolve on different time scales. This result, as 
well as the diversity in mid-IR spectra observed for our 
Class II objects, indicates that age does not seem to be 
the dominating factor determining disk evolution, but 
rather the initial conditions of the star-forming core, like 
its mass and angular momentum. 

While there likely is an evolutionary sequence from a 
classical T Tauri star with an accretion disk to a weak- 
lined T Tauri star with a passive disk to a T Tauri star 
with little or no disk left, there does not seem to be 
an absolute, evolutionary timescale regulating when the 
transitions from one stage to the next occur. Therefore, 
we can witness circumstellar disk evolution in young star- 
forming regions such as Taurus and infer the major steps 
in the process by analyzing a large sample of pre-main- 
sequence stars; this paper is a contribution to this effort 
by presenting and analyzing unprecedented data from 
the Spitzer Space Telescope. 
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Fig. 1. — SED plots of the Class II objects in our sample, ordered alphabetically by their name. In addition to the IRS spectrum, we 
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Fig. 2. — SED plots of the Class III objects i n our sa mple, ordered alphabetically by their name. All fluxes, including the IRS spectrum, 
were dereddened using Mathis's reddening law fMathis 1990) and the extinctions listed in table|3] See Fiaure[J\ legend for a description of 
the plotting symbols. 
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Fig. 10. — Strength of the continuum-subtracted 10-fim feature, normalized to the continuum, vs. the spectral index between 6 and 25 
fim, 116—25- The data points belonging to the groups defined in our morphological sequence are identified by different plotting symbols: 
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Fig. 15. — Median SED from 1.25 to 34 /im computed using 2MASS photometry and the IRS spectrum of only those Class II ob jects with 
spectral types between K5 and M2 (red), as well as the median SED constructed using IRAC data from Hartmann ct al. 1 2005a) (purple), 
compared to two models with M = 3 X 10 — 8 Mq yr — 1 , i=60°, dust consisting of amorphous pyroxene and graphite, and two values of 
settling parameter t (upper, light blue line, e = 0.1; lower, green line, e = 0.001). The dotted line represents the photosphere (SED of the 
WTTS HBC 427). 



T Tauri stars in Taurus Observed with Spitzer IRS 



45 



TABLE 1 

Properties of Observed Class II Objects in Multiple Systems 



Name 


N 


Separation ["] 


Spectral 


CTTS (C)/ 


Brightness 


Ref. 








Type 


WTTS (W) 


Comparison 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


rMQ'yn i oc;crn ( A d\ 
U4o/U+zooy (A,rS) 


2 


yt Off 


t/q ATI 

~ IVo— Ml 




A>B 


20,21 




Z 


i q n" 


MU 




A 3> a 


99 9 

zz,z 


Oolvu lau/o (A,!}^ 


z 


z.UO 


Ml 




A > r> 


1 
1 


P7 Tan f A ~R\ a 


2 


u.oz 


M3 


vv 

v V 


A ~ B 


\ 


V)T) Tan (A T$) a 


2 


0.56" 


ivi o v i o 




A ~ B 


\ 


DF Tan (A R1 


2 


0,09" 






A ~ B 


\ 


T~\TLT T i 

UH lau 


9 


10 


A/f 9j_AAO 
JV1 Z -p 1V1 Z 


ri 
v , 


nn T~1T 

un rjr 


t .o 


Fil-f Tan f A "R^ 
JJJ\ ±aU y±\^L)) 


o 
Z 


Z.OU 


T^Q_l_A/T1 
ivy-|-ivii 




A ^ R 


1 7 


DO Tan (A 


2 


SB 


K5 


Q 


A ~ B 


12 


F041 924-2647 ('A R1 


2 


23. 3" 






A > B 


23 


F04297+2246 (A,B) 


2 


6.6" 






A > B 


23 


FM Tau c 


2 


37.3" 


M2 


c 


FM < V773 


7.8 


FO Tau (A,B) 


2 


0.15" 


M2+M2 


C 


A ~ B 


1 


FQ Tau (A,B) 


2 


0.76" 


M3+M3.5 


c 


A~B 


1,17 


FS Tau (Aa,Ab,B) 


3 


0.23" (Aa,Ab), 20" (A-B) 


M1+M4 (Aa,Ab) 


c 


Aa > Ab 


1,24 


FV Tau (A,B) d 


2 


0.72" 


K5+K6 


c 


A~B, FV>FV/c 


1 


FX Tau (A,B) 


2 


0.89" 


M1+M4 


c+w 


A > B 


1.6 


FZ Tau e 


2 


16.9" 


M0+K5 


c 


FZ > FY 


7.8 


GG Tau (A,B)f 


1 


10.3" 




c 


A>B 


1 


GG Tau (Aa,Ab) 


2 


0.25" 


K7+M0.5 


c 


Aa> Ab 


1 


GG Tau (Ba,Bb) 


2 


1.48" 


M5.5+M7.5 


c 


Ba>Bb 


1 


GH Tau (A,B) 9 


2 


0.31" 


M1.5+M2 


c 


A~B 


1 


GI Tau h 


2 


12.9" 


K6 


c 


GI ~ GK 


5.6 


Tan ( A Rl' 1 


o 
Z 


z.o 


xV ( 


p 


A ^ R 


o.o 


PW Tan ( A R1 


9 
Z 


U.OO 


A/T9 ^ 




A 

. \ 1 > 


1 9fi 

r ,zd 


Ham R ^7 fAa Ah R1 




9 ft?"? A Aa Ah"! 

_.Uz_- 1/1,1) 1, U.OO I jlcl,ii U J 


K7+M1 


Q 


A fl n Ah A^R 


1 11 


RT<" Tan A R1 


9 
Z 




A/T1 _i_A/T9 


p 


A ^ R 
a r_> 


1 7 


HW Tan (A R1 

111N _La,Ll I il - JLj 1 


2 


3.11" 


K"^-UA/T4 

IidTlvlt 




A "S> R 

- i ::' -D 




rrp Ten f A R^* 
il r lau (^-rt.j.Dj 


9 
Z 


U.Ul ( 


W 9 
iVZ 


p 


A \ "R 
a r> 


7 o i c 


TQ TV / A T3\ 

lb lau (A,rSJ 


2 


0.22" 


IV / +1V14.0 


/^i 1 TUT 


A > U 


1 


11 lau (A,dJ 


2 


2.39" 


lvo+M4 






1.6 


RW Aur (A,B,C) J 


3 


1.42 (A-BC), 0.12 (B-C) 


K1+K5 (A,B) 





A>B>C 


1,10 


rr rji /at rt Or.\ 

1 lau (lN,ba,bb) 


3 


U.7U (JN-b), U.l (ba-bb) 


K0 


c 


N ~ Sa ~ Sb 


1,3 


T TV T^i. f A T3 P^ 

Ua lau yJ^^D^j 


1 


efi"^A I5\ 9 ^q"^a p^ 
o.oD (A-rJJ, z.oo (A-L^j 


rv i3 + 1VI z + 1V1 o 


P 1 ^ST" 1 W7 
0-+-W-T-W 


A \ "D A P 

A > rS, A ^> U 


i 
i 


UX Tau (Ba,Bb) 


2 


0.138" 


M2 


w 


Ba> Bb 


11 


UY Aur (A,B) 


2 


0.88" 


M0+M2.5 


c 


A >B 


1,17 


UZ Tau (A,Ba,Bb) 


4 


SB (A), 3.54" (A-Ba), 0.37"(Ba-Bb) 


M1+M2+M2 


c 


A>B, Ba~Bb 


1,13 


V710 Tau (A,B) 


2 


3.17" 


M0.5+M2 


c+w 


A ~ B 


1 


V773 Tau (AB,C,D) 


4 


SB (AB), 0.12" (AB-C), 0.24" (AB-D) 


K2+M0 (AB,C) 


w+c 


D > C> AB 


1,4 


V807 Tau (A,Ba,Bb)f 


3 


0.30" (A-B), 0.04"(Ba-Bb) 


K7+M3 


c+w 


A > B, Ba~Bb 


25,1 


V892 Tau (Aa,Ab,B) 


3 


0.06", 4.10" 


B9+M2 


w 


Aa~Ab, A>B 


16,19 


V955 Tau (A,B) 


2 


0.33" 


K5+M1 


c 


A > B 


1 


VY Tau (A,B) 


2 


0.66" 


M0 


w 


A > B 


1 


XZ Tau (A,B) 


2 


0.30" 


M3+M1.5 


c 


B > A 


1 


ZZ Tau IRS fe 


2 


35" 


M4.5 


c 


ZZ IRS > ZZ 


2 


ZZ Tau (A,B) 


2 


0.04" 


M3 


c 


A >B 


9 



References. (1) IWhite fc Ghed 1200111: (2 ) IWhite fc Hillenbrandl 1200^1: (3) IKoreskd i2000D; ( 41 iDuchene et alj 120031): (5 ) 

IKenvon fc Hartmanri| 119931 : ( 61 IDuchene et all [l99flH: f7~IIMonin e t alj 119931: (81 IHa rtigan. Strom, St roml Il994l> : (9) |^imon et al.| |l990D: 

(10l | Ghez^t^nTTT9*9all: (11 1 IDucheneHl999l1: ( 12l |Mathieu e t*afTTTo97^ : (131 IMathieu. M artlii^^agazzu ni996ll: (14 1 |^Mt^*| r!imimbra^ 

j200ai: (15) I Richichi et alj I1994H: (16) [Leinert m cMchirfc Haad il99 j|l: ( 1 rT^ Mfa^T ^^erlvoTWKm : (18) IPadgett eHuTTT99H|l : (191 
ISmith et ail 120051); (20 1 lltoh e t al.| |199 9fl : (21) lltoliTt^nT2002ll: (221 IDuchene eT^nTT2004ll : (23) 2MASS All-Sky Catalog of Point Sources; 

(241 IMundt et alj 1198411 : (25) IsThaeferTtaTt I2003H : (26) IWhite fc Basril 1200311 

Note. — Column (1) gives the name of the object, column (2) the number of components, column (3) their separation in arcscconds (SB 

means "spectroscopic binary"), column (4) the spectral type (if only one type is listed, it is usually the one of the brighter star in the system), 

column (5) lists whether the components are classical or weak-lined T Tauri stars (if only one identifier is given, it applies to all components), 

column (6) shows how the components compare in brightness, and column (7) gives the references for the data listed in the previous columns. 

a CZ Tau (A,B) is 30" 1 from DD Tau (A,B). 
The other component in the binary is DI Tau, which is itself a 0"l2 binary. 

c The other component in the binary (if it is a real bound system) is V773 Tau, a quadruple system. 

d FV Tau/c (a 0"7 binary) is 12"3 from FV Tau. 

e The other component in the binary is FY Tau. 

f We only observed GG Tau Aa+Ab with the IRS. 

9 GH Tau and V807 Tau are separated by 21" 5. 

h GI Tau and GK Tau are separated by 12"9. 

1 HP Tau might form a triple systems with HP Tau/Gg (2l" 4 away) and HP Tau/G3 (17"3 away). 
J RW Aur C could be a false detection IWhite fc Ghezl200]D . 
k ZZ Tau IRS is 35" from ZZ Tau. 
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TABLE 2 

Properties of Observed Class II Objects 



Name 


Adopted 


A v 


CTTS (C)/ 


Rcf. 




Spectral Type 




WTTS (W) 




(1) 


(2) 


(3) 


(4) 


(5) 


04108+2910 


M0 


1.40° 


C 


1,3 


04187+1927 


M0 




C 


1,4 


04200+2759 






C 


1,5 


04216+2603 


Ml 




C 


1,3 


04303+2240 




11.7 


C 


6 


04370+2559 (A,B) 




9.82 


C 


7 


04385+2550 (A,B) 


M0 


7.80 a 


C 


6,7 


AA Tau 


K7 


1.75 a 


C 


1,3 


AB Aur 


AO 


0.25 


C 


8,3 


BP Tau 


K7 


1.00" 


C 


1,3 


CI Tau 


K7 


2.00 a 


C 


1,3 


CoKu Tau/3 (A,B) 


Ml 


5.00 a 


C 


1,2 


CoKu Tau/4 


Ml 


3.0 


w 


9,3 


CW Tau 


K3 


2.75 a 


C 


1,3 


CX Tau 


M0 


1.30" 


C 


1,3 


CY Tau 


K7 


1.70 a 


C 


1,3 


CZ Tau (A,B) 


Ml 


2.45 a 


w 


1,3 


DD Tau (A,B) 


M3 


1.00" 


C 


1,2 


DE Tau 


M0 


1.20° 


C 


1,3 


DF Tau (A,B) 


M0 


1.60° 


C 


1,2 


DC Tau 


K6 


1.60 


C 


1,10,3 


DH Tau 


M0 


1.70 a 


C 


1,3 


DK Tau (A,B) 


M0 


1.30 a 


C 


1,2 


DL Tau 


K7 


1.50° 


c 


1,2 


DM Tau 


Ml 


0.72 a 


c 


1,2 


DN Tau 


M0 


0.60 a 


c 


1,3 


DO Tau 


M0 


2.05 a 


c 


1,3 


DP Tau 


M0 


0.60 a 


c 


1,3 


DQ Tau (A,B) 


M0 


1.60 a 


c 


1,3 


DR Tau 


K7 


1.20 


c 


1,3 


DS Tau 


K5 


1.10 a 


c 


1,3 


F04101+3103 


Al 


1.90" 


c 


1,4 


F04147+2822 


M4 


2.50 a 


c 


1,3 


F04192+2647 A 










F04262+2654 










F04297+2246 A 










F04570+2520 










FM Tau 


M0 


1.40 a 


c 


1,3 


FN Tau 


M5 


1.35 


c 


1,3 


FO Tau (A,B) 


M2 


3.03 


c 


2 


FP Tau 


M4 


0.00 


c 


2 


FQ Tau (A,B) 


M2 


1.87 


c 


1,3 


FS Tau (Aa,Ab) 


Ml 


1.43 a 


c 


1,2 


FT Tau 


C 






1 


FV Tau (A,B) 


K5 


5.33 


c 


2 


FX Tau (A,B) 


Ml 


2.00 a 


c 


1,3 


FZ Tau 


M0 


3.70 a 


c 


1,3 


GG Tau (Aa,Ab) 


M0 


1.00 a 


c 


1,2 


GH Tau (A,B) 


M2 


0.97 a 


c 


1,2 


GI Tau 


K6 


2.30 a 


c 


1,3 


GK Tau (A,B) 


M0 


1.15 a 


c 


1,3 


GM Aur 


K3 


1.21 


c 


2 


GN Tau (A,B) 


M2 


3.50 a 


c 


11,12 


GO Tau 


M0 


2.02 a 


c 


1,3 


Haro 6-13 


M0 


11.9 


c 


6 


Haro 6-37 (Aa,Ab,B) 


K7 


3.80 a 


c 


1,2 


HK Tau (A,B) 


Ml 


2.70 a 


c 


6 


HN Tau (A,B) 


K5 


1.50 a 


c 


1,2 


HO Tau 


M0 


1.30° 


c 


1,3 


HP Tau (A,B) 


K3 


2.80 a 


c 


1,3 


HQ Tau 










IP Tau 


M0 


0.51 


c 


2 


IQ Tau 


Ml 


1.44 


w 


2 


IS Tau (A,B) 


M0 


3.22 a 


c 


1,2 


IT Tau (A,B) 


K2 


3.80 a 


c 


1 


LkCa 15 


K5 


1.20 a 


c 


1,3 


MHO 3 


K7 


8.30 a 


c 


1,13 


RW Aur (A,B,C) 


K3 


0.50 a 


c 


1,2 


RY Tau 


Gl 


2.20 


c 


14,3 


SU Aur 


Gl 


0.90 


c 


14 


T Tau (N,Sa,Sb) 


K0 


1.75 a 


c 


1,2 


UX Tau (A,Ba,Bb,C) 


K5 


0.70 a 


c 


2 


UY Aur (A,B) 


K7 


2.10 a 


c 


1,3 
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TABLE 2 — Continued 




- \ V 1 1 ' j / Lv A 1 


J v V 


CTTS (d\ 1 


Ref. 




Slr^pr'tva 1 Twc^o 




WTTS f"Wl 






(1) 




V *J 




UZ Tau (A,Ba,Bb) 


Ml 


1.00 


c 


2 


V410 Anon 13 


M6 


5.80 


c 


1,15,2 


V710 Tau (A,B) 


Ml 


1.90 a 


c 


1,2 


V773 Tau (AB,C,D) 


K3 


2.00 a 


w 


1.2 


V807 Tau (A,Ba,Bb) 


K7 


0.60 a 


c 


1.2 


V836 Tau 


K7 


1.10 a 


c/w 


1,3 


V892 Tau (Aa,Ab,B) 


B9 


8.00 a 


w 


16,3 


V955 Tau (A,B) 


K5 


3.72 


c 


2 


VY Tau (A,B) 


MO 


1.35 a 


w 


1,3 


XZ Tau (A,B) 


M2 


2.90 a 


c 


1.2 


ZZ Tau IRS 


M5 


1.50 a 


c 


1,6 


ZZ Tau (A,B) 


M3 


1.44 a 


c 


1,3 



Refe renc es. — (l) IKenvon fo Hartmanr] 1199511: (2HWhite fe Ghe3 120011): OIIKenvon et all 119931: (4irKenvon et alj 119901): C5~) IKenvon et*afl 
I1994a[): f6)IWhite fcHmenbrandl 120041): (71 | Itoh et alj 120021): (8)IDeWarf et al.1 1200311: (9) ID'Alessio et alj I2005|): aollGullbring et alj 120001): 
fllllwhite fc BasriH2003D : ( 12) ILuhmanH2004h : fl3) IBriceno et al.Nl99gf ): (14) ICalvet et all 120041) : f!5) IFurlan etall I2005aj) : (16) ISmith et alj 
12003) 

Note. — Column (1) gives the name of the object, column (2) the adopted spectral type, column (3) the adopted visual extinction (see also 
note (a) below), column (4) lists whether the object (the primary in case of multiple systems) is a classical or weak-lined T Tauri star, and column 

(5) gives the references for the data listed in the previous columns. 

a Valu e of Ay derived from t he observed VT colors, spectral types, and intrinsic colors for main-sequence stars listed in Kcnvon & Hartmann 
using the lMathisI I199CD reddening law. In some cases the derived values of Ay were adjusted to yield a good match between the dereddened 
and true photospheric colors. 
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TABLE 3 

Properties of Observed Class III Objects 



Name 


Multiplicity" 


Adopted 
Spectral Type 


Ay 


CTTS (C)/ 
WTTS (W) 


Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Anon 1 


s 


]\ TO 


Z. / U 


V V 


1,0 


JJ1 lau 


U.lz , lb 


A in 
JVLU 


0.83 b 


VV 


7,2,5 


r r lau i\ , J3 ^ 


U.UZO 


k"7 
i\ t 


z.zz 


W7 
VV 


O A 


HBC 35b 


1.33" 


K3 


0.42 b 


W 


3,2 


TJ 73 /~> goo 

not ooo 


s 


i.- r i 
ivl 


0.0 


\ K r 
VV 


1 


TTT3 /~> QOO 

HBC oyz 


s 


1/ r 


0.05 & 


\ K t 
VV 


2,5 


IlDb 4Zo 


U..Z4 , 1U.4 


TV/T1 

M 1 


3.40 b 


VV 


1,2 


TJT2 /~* yir|7 

HBC 4Z/ 


bB [rsj U.Uo j 


I." rt 

Kb 


0.50 b 


A A ; 

VV 


2,5,9 


xlJJ ZOOO / Z 


s 


bo 


r\ AAb 
U.44 


VV 


2,5 


TJD TV,,, /Oof 

nr lau/CrZ J 


y.y 


CO 


2.63 & 


VV 


2,4,5 


Hubble 4 


s 




2.40^ 


VV 


1 


1W lau (A, 15 J 


u.zy 




i nab 
l.Ub 


VV 


i 
1 




S 


K7 


n n 


V V 


1,0 


LkCa 1 


ft 


M4 


0.0 


VV 


1 


LkCa 3 (A,B) 


0.48" 


Ml 


0.88 & 


vv 


2 1 


LkCa 4 


s 


K7 


1.10 b 


VV 


1 


LkCa 5 


s 


M2 


0.40 & 


w 


1 


LkCa 7 (A,B) 


1.02" 


K7 


1.00 6 


w 


1 


LkCa 21 


s 


M3 


0.30 b 


VV 


2,6 


V410 Tau (A,B,C)s 


0.07", 0.29" 


K3 


0.94 & 


VV 


1 


V410 X-ray 3 


s 


M6 


0.80 & 


VV 


1 


V819 Tau 


10.5" 


K7 


1.70 6 


w 


1,8 


V826 Tau 


SB 


K7 


0.70 b 


w 


8,5 


V827 Tau 


s 


K7 


1.00 6 


VV 


1 


V830 Tau 


s 


K7 


0.50 b 


VV 


1 


V928 Tau (A,B) 


0.19" 


M0.5 


2.10 & 


w 


2,6 



Ref ere nces. — (l)lWhite fc Ghed I200U) (2)lKenvon fc Hartmannl il993) OI lDuchene et all Il99s| f4l lSimon et alJ 119951) l^ lKenvon et all 
1199811 f6^ lMartin et all i!994l) m iGhez eHifl 11995ft C8i lLeinert et all 11993ft O USteffen et all i20QlD 

Note. — Column (1) gives the name of the object, column (2) the multiplicity of the object (see footnote a), column (3) the adopted spectral 
type, column (4) the adopted visual extinction (see also note (b) below), column (5) lists whether the object (the primary in case of multiple 
systems) is a classical or weak-lined T Tauri star, and column (6) gives the references for the data listed in the previous columns. 
a The letter "s" means single star; for multiple systems, the separation between the components in arcscconds is listed; SB means "spectroscopic 

binary" . 

h Valu e of Ay- d erived from th e observed V-I colors, spectral types, and intrinsic colors for main-sequence stars listed in Kcnvon Sc Hartmann 
11993) . using t.he lMathil 1199(11 reddening law. In some cases the derived values of Av were adjusted to yield a good match between the dereddened 
and true photospheric colors. 

c DI Tau is separated by 15" from DH Tau; DI Tau is a 0"l2 binary. 
d The other component is HBC 357. 

e HBC 423 is also knows as LkHa 332/Gl. The separation between components A and B of HBC 423 is 0. 24; the separation between HBC 423 
and V955 Tau is 10'/4. 

f HP Tau/G2 is likely part of a triple system consisting of HP Tau, HP Tau/G2, and HP Tau/G3; the separation listed is the distance between 
HP Tau/G2 and HP Tau/G3 (also a WTTS). 

9 The value of 0"07 is the separation between A and B, 0"29 between A and C. 
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TABLE 4 

Spectral Index ne-25 and 10-fiu feature strength for the Class II 
Objects in our Sample 



Name 


7^6 — 25 


I _T 1(1 r mrt ) / _T i-i-.n + 

1U COilL// ■*- COIll 


Name 




t r in F f-nnt ) / .T rrin+ 

_LU COilL// - 1 - COilL 


04108+2910 


-0.90 


0.01 


FS Tau 


-0.15 


0.26 


04187+1927 


-0.56 


0.28 


FT Tau 


-0.48 


0.43 


04200+2759 


-0.32 


0.44 


FV Tau 


-0.62 


0.17 


04216+2603 


-0.57 


0.27 


FX Tau 


-0.58 


0.81 


04303+2240 


-0.70 


0.21 


FZ Tau 


-0.86 


0.18 


04370+2559 


-0.24 


0.61 


GG Tau 


-0.47 


0.64 


04385+2550 


0.35 


0.31 


GH Tau 


-0.75 


0.26 


AA Tau 


-0.65 


0.35 


GI Tau 


-0.47 


0.51 


AB Aur 


-0.05 


1.08 


GK Tau 


-0.20 


1.03 


BP Tau 


-0.48 


0.55 


GM Aur 


0.46 


1.19 


CI Tau 


-0.54 


0.54 


GN Tau 


-0.75 


0.59 


CoKu Tau/3 


-0.95 


0.41 


GO Tau 


-0.37 


0.40 


CoKu Tau/4 


0.95 


0.85 


Haro 6-13 


0.15 


0.35 


CW Tau 


-0.83 


0.16 


Haro 6-37 


-0.87 


0.18 


CX Tau 


-0.14 


0.48 


HK Tau 


0.09 


0.39 


CY Tau 


-1.22 


0.12 


HN Tau 


-0.45 


0.51 


CZ Tau 


-0.01 


0.53 


HO Tau 


-0.66 


0.58 


DD Tau 


-0.64 


0.21 


HP Tau 


-0.07 


0.48 


DE Tau 


-0.44 


0.43 


HQ Tau 


-0.50 


0.70 


DF Tau 


-1.13 


0.17 


IP Tau 


-0.45 


0.99 


DC Tau 


0.00 


0.00 


IQ Tau 


-1.12 


0.35 


DH Tau 


-0.67 


0.70 


IS Tau 


-0.85 


0.44 


DK Tau 


-0.72 


0.83 


IT Tau 


-1.05 


0.14 


DL Tau 


-0.67 


0.09 


LkCa 15 


-0.36 


1.47 


DM Tau 


0.72 


0.96 


MHO 3 


0.94 


0.69 


DN Tau 


-0.45 


0.22 


RW Aur 


-0.50 


0.32 


DO Tau 


-0.27 


0.18 


RY Tau 


0.05 


1.36 


DP Tau 


-0.23 


0.35 


SU Aur 


0.37 


0.98 


DQ Tau 


-0.38 


0.14 


T Tau 


-0.04 


-0.10 


DR Tau 


-0.61 


0.27 


UX Tau A 


-0.08 


0.15 


DS Tau 


-0.90 


0.55 


UY Aur 


0.08 


0.40 


F04101+3103 


0.75 


1.03 


UZ Tau 


-0.66 


0.45 


F04147+2822 


-0.75 


0.56 


V410 Anonl3 


-0.61 


0.20 


F04192+2647 


-0.67 


0.28 


V710 Tau 


-0.87 


0.21 


F04262+2654 


-0.38 


0.19 


V773 Tau 


-0.76 


0.09 


F04297+2246 A 


-0.24 


0.98 


V807 Tau 


-0.91 


0.16 


F04570+2520 


-1.04 


0.25 


V836 Tau 


-0.86 


0.79 


FM Tau 


-0.19 


0.72 


V892 Tau 


1.28 


0.65 


FN Tau 


-0.04 


0.39 


V955 Tau 


-0.78 


0.24 


FO Tau 


-0.52 


0.26 


VY Tau 


-0.63 


0.49 


FP Tau 


-0.39 


0.22 


XZ Tau 


-0.23 


0.09 


FQ Tau 


-0.78 


0.19 


ZZ Tau IRS 


-0.02 


0.31 








ZZ Tau 


-0.94 


0.26 



Note. Spectral indices are defined as n — d\og(\F\) f d\og(\) . Sec text for the definition of the strength of the 10-/i.m feature, 

(F 10 — F cont ) /F cont . 



50 



Furlan et al. 



TABLE 5 

Spectral Indices of the Class II and III Objects in our Sample 



Name 




Tl\^ 25 


Name 


7lQ — 13 


71,13 25 


Class II Objects 












04108+2910 


-1.17 


-0.55 


GN Tau 


-0.89 


-0.63 


04187+1927 


-0.49 


-0.66 


GO Tau 


-0.80 


0.17 


04200+2759 


-0.47 


-0.11 


Haro 6-13 


-0.56 


0.90 


04216+2603 


-0.75 


-0.33 


Haro 6-37 


-1.17 


-0.53 


04303+2240 


-0.80 


-0.74 


HK Tau 


-0.72 


1.12 


04370+2559 


-0.37 


-0.22 


HN Tau 


-0.70 


-0.14 


04385+2550 


0.13 


0.53 


HO Tau 


-0.89 


-0.37 


AA Tau 


-0.87 


-0.37 


HP Tau 


-0.43 


0.37 


AB Aur 


-0.79 


0.93 


HQ Tau 


-0.76 


-0.17 


BP Tau 


-0.67 


-0.24 


IP Tau 


-0.94 


0.20 


CI Tau 


-0.95 


-0.02 


IQ Tau 


-1.30 


-0.90 


CoKu Tau/3 


-1.04 


-0.91 


IS Tau 


-0.87 


-0.89 


CoKu Tau/4 


-0.60 


2.98 


IT Tau 


-1.38 


-0.67 


CW Tau 


-1.26 


-0.30 


LkCa 15 


-1.30 


0.87 


CX Tau 


-0.37 


0.14 


MHO 3 


0.72 


1.10 


CY Tau 


-1.45 


-0.95 


RW Aur 


-0.67 


-0.27 


CZ Tau 


0.42 


-0.61 


RY Tau 


-0.16 


0.30 


DD Tau 


-0.73 


-0.53 


SU Aur 


-0.32 


1.26 


DE Tau 


-0.87 


0.12 


T Tau 


-0.56 


0.63 


DF Tau 


-1.36 


-0.85 


UX Tau A 


-1.99 


2.45 


DC Tau 


-0.20 


0.24 


UY Aur 


-0.01 


0.17 


DH Tau 


-1.69 


0.66 


UZ Tau 


-0.92 


-0.34 


DK Tau 


-1.00 


-0.37 


V410 Anon 13 


-0.91 


-0.31 


DL Tau 


-0.78 


-0.54 


V710 Tau 


-1.07 


-0.63 


DM Tau 


-0.11 


1.82 


V773 Tau 


-1.10 


-0.34 


DN Tau 


-0.60 


-0.27 


V807 Tau 


-1.61 


0.01 


DO Tau 


-0.56 


0.08 


V836 Tau 


-1.37 


-0.20 


DP Tau 


-0.43 


0.02 


V892 Tau 


1.75 


0.53 


DQ Tau 


-0.47 


-0.28 


V955 Tau 


-0.89 


-0.68 


DR Tau 


-0.88 


-0.27 


VY Tau 


-1.11 


-0.02 


DS Tau 


-1.03 


-0.73 


XZ Tau 


-0.35 


-0.10 


F04101+3103 


0.97 


0.43 


ZZ Tau IRS 


-0.42 


0.49 


F04147+2822 


-0.83 


-0.70 


ZZ Tau 


-1.15 


-0.68 


F04192+2647 


-1.10 


-0.10 








F04262+2654 


-0.89 


0.30 


Class III Objects 






F04297+2246 A 


-0.41 


-0.01 


Anon 1 


-2.78 


-2.71 


F04570+2520 


-0.79 


-1.37 


FF Tau 


-2.76 


-2.70 


FM Tau 


-0.45 


0.14 


HBC 388 


-2.78 


-3.02 


FN Tau 


-0.29 


0.27 


HBC 427 


-2.85 


-0.70 


FO Tau 


-0.82 


-0.18 


HD 283572 


-2.85 


-2.95 


FP Tau 


-0.84 


0.21 


IW Tau 


-2.82 


-2.69 


FQ Tau 


-1.12 


-0.36 


LkCa 3 


-2.81 


-2.65 


FS Tau 


-0.56 


0.37 


LkCa 4 


-2.74 


-2.72 


FT Tau 


-0.62 


-0.31 


LkCa 5 


-2.77 


-3.06 


FV Tau 


-0.83 


-0.42 


LkCa 7 


-2.94 


-2.46 


FX Tau 


-1.05 


0.02 


V410 Tau 


-2.78 


-5.19 


FZ Tau 


-1.14 


-0.56 


V410 Xray 3 


-2.89 


-1.61 


GG Tau 


-0.79 


-0.05 


V819 Tau 


-2.57 


-0.67 


GH Tau 


-1.26 


-0.10 


V826 Tau 


-2.80 


-2.68 


GI Tau 


-0.55 


-0.40 


V827 Tau 


-2.77 


-3.13 


GK Tau 


-0.45 


0.12 


V830 Tau 


-2.65 


-2.99 


GM Aur 


-0.91 


2.28 


V928 Tau 


-2.76 


-2.72 



Note. — Spectral indices are defined as n — d log(AF A )jd log(A). 



T Tauri stars in Taurus Observed with Spitzer IRS 
TABLE 6 

IRS Taurus Median (K5-M2 spectral types) 



Wavelength 


Median 


Lower Quartilc 


Upper Quartilc 


H 


[\og(uF v )] a 


{\og(vF„)] a 


[\og{uF 1/ )] a 


l.ZD 


o mo 
-y.Uoy 


-y.uyz 


-y.uio 


1.65 6 


-9.035 


-9.035 


-9.035 


2.17 6 


-9.163 


-9.195 


-9.114 


5.70 


-9.764 


-9.879 


-9.614 


7.10 


-9.873 


-10.060 


-9.735 


8.00 


-9.905 


-10.081 


-9.712 


9.20 


-9.789 


-9.967 


-9.552 


9.80 


-9.793 


-9.942 


-9.519 


11.30 


-9.902 


-10.060 


-9.615 


12.30 


-10.034 


-10.234 


-9.763 


13.25 


-10.060 


-10.304 


-9.787 


16.25 


-10.066 


-10.237 


-9.777 


18.00 


-10.030 


-10.182 


-9.725 


21.00 


-10.037 


-10.204 


-9.722 


25.00 


-10.093 


-10.278 


-9.778 


30.00 


-10.185 


-10.404 


-9.839 


34.00 


-10.229 


-10.461 


-9.865 



Note. — a Units of median, lower and upper quartiles arc ergs s 1 cm 2 . 
The median at these 3 wavelengths was constructed from reddening-corrected 2MASS data. 



